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This thesis is a three-dimensional numerical analysis of Gas Turbine Test Cell #1 
at Lemoore Naval Air Station, California. The Parabolic, Hyperbolic, or Elliptic 
Numerical Integration Code (PHOENICS) is used to determine the steady-state 
aerothermal characteristics within the test cell during the full power run-up of a GE 
F404 gas turbine engine (afterburner in operation). The method in which PHOENICS 
arrives at a solution is discussed as well as how the code was applied to this problem. 

The solution of the problem is presented in tabular form and the results are 


discussed. Recommendations for the future application of PHOENICS to this picjcet 
are given. 
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I. INTRODUCTION 


Critical in the rework and overhaul of aviation gas turbine engines is the testing 
of the powerplants prior to their return to flight status. This testing consists of “run- 
up’ and engine parameter checks either in the airframe or attached to a test stand. A 
wide variety of test facility designs are currently used by the Navy/ Marine Corps repair 
activiues. Tliese test facilities vary from open-air run-up pads where the aircraft is tied 
down during testing to more complex structures such as hush houses, which allow fully 
enclosed in-airframe testing and fully enclosed test cells where engines are mounted on 
test stands during run-up. The choice of the actual design installation is driven by 
several factors - among them being the number of operating personnel available, the 
type of powerplants to be tested, noise and pollution standards to be met, and cost and 
size constraints. 

Until now, gas turbine test facilities have been designed and built primarily on 
the basis of the intuition and experience of design engineers rather than with scientific 
modeling and testing of proposed plans. This practice resulted in a myriad of different 
structures without necessarily producing the best designs. Most of these test facilities 
performed adequately until the emergence of the new, high powered aviation gas 
turbine engines. The tremendous thrust generated by these engines resulted in 
velocities, pressures, temperatures, noise and air pollution levels too great for the 
current test cells to control [Ref: 1]. A heightened public awareness of pollution and 
noise hazards has led to pressure to closely examine and inprove test sites. The Navy 
Environmental Health Center reported hearing loss claims due to excessive noise at 
$114 million through FY81, with a steady 9% increase annually. Also in 1981, the 
Federal Court decided that military jet test sites must comply with state and local 
resmarions pursuant to Section 233 othe Clean Air Act. [Ref. 2: p. 3] 

Currently there are 75 Navy/Marine Corps activities operating 76 enclosed test 
cells and 5 hush houses. Of the test cells, less than 80% are fully operational and two- 
thirds of them are at least 25 vears old. An analysis conducted in F Y83 concluded that 
costs to rehabilitate these facilities at $84 million. Considering the climbing costs of 
maintainance and repair of the cells, the lack of standardized designs, the safety 


hazards posed, and the pollution’/noise problems, the Department of the Navy has 


chosen to undertake a comprehensive program of research, development, testing and 
evaluation (RDT&E) leading to the construction of standardized test facilities at 
Navy/ Marine Corps aviation repair activities worldwide. [Ref 2: p. 3] 

A portion of this research is the detailed analysis of the existing gas turbine test 
facilities. Information is needed as to the aerothermal conditions within the structure as 
well as the ability of the structure to contain noise and minimize pollutants. In this 
case the PHOENICS code, a FORTRAN-type computer program, is used to analyze 
the performance of Lemoore Naval Air Station Test Cell #1 during the run-up of a GE 
F404 gas turbine engine. A variety of physical parameters, including pressure, 
temperature, velocity, density, enthalpy and direction, are found using PHOENICS. 
This data will then be used by engineers performing the design phase of the RDT&E to 


identify the positive and negative features of the existing facilities. 
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Il. PROBLEM FORMULATION 


A. DESCRIPTION OF THE TEST CELL 

The test cell being investigated in this thesis is formally called the 
Coanda/ Refraction Noise Suppression System. It is one of two of this type built at 
Lemoore Naval Air Station, California. The cells are used in support of the overhaul 
of aviation gas turbines at the depot level. One of the most common engines tested at 
Lemoore and the one used in this project, is the GE F404. The F404 gas turbine is the 
powerplant for the Navy F/A-I8. 

In the normal maintainance cycle of the F/A-18 and when serious engine 
problems arise, the engines are given a major overhaul. The completion of the overhaul 
is marked by the successful run-up of the engine in a test cell such as the Coanda. 
These test cells are subject to heavy use, typically with useage in the order of several 
times weekly and daily when work loads are high. The test of the F404 engine lasts 
approximately 45 minutes, of which the final 15 minutes are at or near full power. 

The Coanda-Refraction Noise Suppression Svstem is a reinforced concrete 
structure which fully encloses the engine test stand. At one end, large double doors 
open to allow the engine to be easily moved in and out of the facility. Once in place, 
tae cneimetis seeuted to the test eradlefthe doors are closed, and the fuel, electrical, 
and instrumentation connections made. The engine is controlled during the test by 
personnel inside the control booth, an acoustically insulated compartment which also 
houses all monitoring and control equipment. In the test room, where the engine 1s 
located, an installed firefighting system as well as peripherals such as observation 
mirrors, hose reels and wash down equipment are found. 

The general dimensions of the Coanda test cell are 30 meters long by 7 meters 
wide. The highest point of the cell, the exhaust stack, is 17 meters. The complex 
geometry of the Coanda cell is shown in Figure 2.1 . A single primary air inlet directs 
air into the cell for ingestion by the engine. This air is divided into two separate 
channels before turning vanes redirect the flow toward the engine intake. Two 
secondary air inlets entrain air for cooling of the engine exhaust gases. The forward 
secondary air inlet allows cooling air to be mixed with the jet blast in the ejector tube, 


a stcel duct encircled with three intake vanes. The after secondary inlet entrains cool 


air Which passes over the Coanda surface and mixes with the flow at the base of the 
exhaust. The mixed primary and secondary flows pass up the exhaust stack to the 
atmosphere. Noise attenuators are placed at the primary and secondary inlets while 
acoustic dampening material lines the walls of the exhaust stack. The concept of the 
svstem is depicted in Figure 2.2 . Screens cover the primary inlet in an attempt to hmut 
the foreign matter ingested into the cell. Use of the cell after idle periods must begin 


with the removal of birds which tend to roost in the inlets. 


B. PROBLEMS AND QUESTIONS ABOUT THE COANDA CELL 

While the Coanda test cells are currently perfornung adequately, there are 
problems associated with them. Investigation of the Coanda Test Cell #1 showed signs 
of moderate spalling and cracking of the concrete walls and there are reports from 
Operating personnel of extreme vibration in the ejector tube and Coanda sur.ace. The 
vibration problem has necessitated the inclusion of structural “bolt-tightening” as part 
of the pre-operation check-off list. 

No specific complaints about excessive noise or pollutants were recorded during 
the on-site inspection, however, the extremely remote locale of the test cells probably 
limits the observation by non-interested personnel. 

There is a lack of information concerning the efficiency of the Coanda test cells. 
While 1t appears that the cell is performing as designed, no phvsical data is available to 
help in confirmation. The object of this thesis was to develop physical data on the 
Coanda cell using purely numerical means. Three-dimensional flow velocities, pressure 
and temperature distribution, along with density and kinetic energy is vital information 
in understanding the test facility and the measure of its performance. 

The approach used to achieve this numerical analysis was through the use of the 
PIEIOENICS computer code. The code, which 1s "discusscd=in more cetan in Charen 
I{I, solved for the unknown physical parameters using only the conditions known to 
exist during the steady state operation of the cell. The information derived from this 
study will be used to determine which features of the Coanda cell mav be useful in 


future test cell designs. 
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lil. THE USE OF PHOENICS IN THE PROBLEM 


A. A BRIEF EXPLANATION OF PHOENICS 
PHOENICS, the computer code used in this analysis of the Coanda test cell, is a 
copyrighted product from CHAM (Concentration Heat and Momentum) Limited. 
Cham is a software and consultancy organization founded in 1972 by D. Brian 
Spalding which specializes in the computer modeling of processes involving fluid flow, 
heat transfer, chemical reaction, and combustion. PHOENICS, which stands for 
Parabolic, Hyperbolic, or Elliptic Numerical Integration Code Series, is a general 
purpose computer code system for simulating the above processes in engineering 
equipment or the environment. [Ref. 3: p. I] 
The basic capabilities of PHOENICS are the simulation of fluid flow and heat 
transfer which are: 
® one, two or three dimensional 
e steady or transient 
e single, two or multi-phase mixtures 
e laminar or turbulent 
e subsonic, supersonic, or transonic 
e chemically reactive or inert 
e compressible or incompressible 
e parabolic, elliptic or hyperbolic 
The user also has command of up to 25 dependent variables, depending on the 
requirements of the problem. Some of these variables include pressure, three 
dimensional velocity components, enthalpy, turbulence, kinetic energy, and radiation 
Were, (Ref 3:.p* 2] 
PHOENICS can be utilized on projects of varying scales and degrees of accuracy. 
The code lends itself to the solution of problems of the “macro” scale as well as those 
of the “micro” scale. The quality of the solution is highly dependent upon the accuracy 
of the initial modeling and the boundary/initial conditions imposed. There 1s a trade- 
off however, a search for more accuracy usually results in more computer tume and 


thus, more expense. 


The problem to be solved is established in the “QI” file. This file, which is created 
entirely by the user, defines all aspects of the problem including physical dimensions, 
boundarv/initial conditions, grid spacing, porosities, and known values of temperature, 
velocity, pressure, etc. The QI data input file is then read and interpreted by the 
“SATELLITE” program into a language recognizable by EARTH, the main body of 
the PHOENIGS code: 

The method in which PHOENICS arrives at the desired quantities is through the 


solution of the following generalized equation which governs conservation law 


phenomena: 
G ee a 
Ge (rp@) + div(rpV@-rl geradg) = Se 
where: 
r = phase volume fraction 
Dp = dependent variable 
lo = exchange coefficient (laminar or turbulent) 
S¢ = source or sink term 
O = density 
V = velocity vector 


= time 


ce 
| 


The users choice of the property that @ will represent determines which conservation 
equation the “generalized equation” will represent (1-e. continuity equation, heat 
transfer equation, momentum equations, etc). [Ref. 3: p. 6] 

The distribution of the values in the grid is computed through a finite-domain 
(FED) method. This method involves the solution of a set of algebraic relationships of 
the following type: 


ap Or ss ay; D yy, te arn Px) is acD¢ a apyPry = ap Py a5 ay P+ oS 
DP = ne 


ap + aw + aQ tT ag + ayy + ay + Ay F AY 


where: 
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@D = the dependent variable in question 


P = the point in the grid in question 
= subscript indicating the cell to the EAST 
W = subscript indicating the cell to the WEST 
N = subscript indicating the cell to the NORTH 
S = subscript indicating the cell to the SOUTH 
H = subscript indicating the cell that is HIGH (above) 
le = subscript indicating the cell that is LOW (below) 
T = subscript indicating the earlier time value (transient only) 


S and ap express the influence of a source of the entity ® while the “a” represents 
coefficients with the dimensions of mass per unit time. The values of all parameters 
except velocity are calculated at cell centers. The three velocity components (U,V, and 
W) are calculated at the cell boundaries. [Ref. 4: p. 2.5] 

The solution of a grid of even a small number of cells is very complicated, with 
numerous, strongly coupled algebraic equations to be solved. PHOENICS uses an 
iterative “guess and correct” process with adjustments and hopefully, improvements, 
made at every iteration. 

These iterations are performed in “sweeps”, through any of a combination of 
possible schemes. These schemes, which are described in detail in the PHHOENICS 
literature, include “slabwise”, “whole-field”, and “point-to-point” solutions. Correct 
application of any of these schemes to a specific problem will result in identical 
solutions, however the difference lies in computer time and storage necded for the 
different methods. It should be noted that certain chosen variables can be found using 
one scheme while other variables may be solved using a different scheme... in the same 
problem. The advantage of using this approach is the possible reduction of computer 
time needed and quicker convergence of the problem to a correct solution. [Ref. 4: p. 
25 

Convergence (or divergence) of the solution is indicated by “residual” values 
given for each unknown at every sweep. These residuals will decrease as the problem 
progresses if convergence is occuring. PHOENICS has arrived at a solution when the 
residuals no longer decrease. The value of this residual may be calculated by the user 
using the formula described in the PIIOENICS literature. [t should be noted, however, 


that the convergence is rarely monotonic due to the strong coupling between the 
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unknowns and the varying speeds at which unknowns will converge. Steadily 
increasing residual values indicate divergence and corrective methods should be 
employed. These methods are also thoroughly discussed in the PHOENICS literature. 

The matter of problem divergencc is one which illustrates that the successful use 
of PHOENICS is highly dependent upon experience and a “feel” for the behavior of the 
code in a given situation. One of the more powerful tools used in controlling 
divergence is the RELAX command. Relaxation placed on a dependent variable 
controls its variation in consecutive sweeps, thereby limiting the possibility of 
divergence to an excessively large or small value. Smaller RELAX values result in 
smallcr variation, with very small values effectively “freezing” the value. At this point it 
becomes apparent that if the RELAX values are too large, divergence will continue 
while if the valucs are too small, convergence will require excessive computer time. 

For problems of moderate size and larger, the EARTH program requires 
considerable computer time to run. The complexity of the computations, and sizc of 
the problem dictate the actual time required. Data from successful runs is stored in a 
“DFO9” file which is extremely useful as a starting point for subsequent runs, or 
“restarts.” Depending upon the definition of the problem, PHOENICS may require 
several hundreds or thousands of sweeps to achieve convergence. 

Typically, the routine followcd using the PHOENICS code was one of the 
submussion of overnight runs in the order of 250-300 sweeps. The rcsults were reviewed 
the following morning and any adjustments to the code were made. Depending on the 
size of the code and the complexity of the problem, a fcw small runs of 50 sweeps 
could be made during the day. The submission of large runs, however, was not 
possible until the ncxt “night-time rate” was available. This routine was followed until 


the convergence of the problem was achieved. 


B. MODELING OF THE PROBLEM USING PHOENICS 

The modeling of the test cell began with an analysis of the physical dimensions of 
the structure. The blueprints were used to accuratcly determine the size and shape of 
the cell as well as the position of the engine test stand and other internal features. On 
site inspection of the facility revealed those aspects, such as clearances and complex 
surface shapes, that were not readily apparcnt from the blueprints. The next step in 
the understanding of the physical geometry of the problein was thc representation of 
the structure on a cartesian coordinate system. The cell was carefully redrawn to scale 


and divided into appropriate subdivisions in the X,Y and Z directions. Since the 
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PHOENICS code uses the SI system of measurement, all lengths and parameters had 
to be converted from the english system. Lengths were converted with an accuracy of 
0.001 meters. Geometric shapes were approximated as accurately as possible, with 
deviations in Situations such as square vice rounded corners and rectangular vice 
elliptical tunnels, etc. The centerline of the test cell was used as a line of symmetry and 
the model was constructed so that only half the actual test cell was modeled. This 
technique allowed a large reduction in the computer code needed for the model. 

The structure was modeled in the PHOENICS code by setting boundaries and 
conditions in the “QL” file. Statements in this file describe the physical dimensions and 
the known boundary conditions. It should be noted that PHOENICS considers any 
Known value within the problem to be a “boundary condition”. The established 
cartesian coordinate system was used to precisely describe the relative positions of the 
walls, deck, openings, and “patches” where particular conditions are known to exist. 
The Z-axis was oriented along the major flow direction, as suggested by CHAM. For 
the model analyzed in this thesis, the number of cells in the X,Y, and Z directions were 
6, 16 and 19, respectively. Grid size was determined by three factors: test cell geometry, 
problem complexity, and program storage dimension available. 

An overlay of the grid onto the cell diagram is shown in Figures 3.1 and 3.2. In 
these figures, it is readily apparent that structures within the cell dictated the absolute 
minimum number of grid spacings required. This is due to the fact that test cell 
boundaries, structures and conditions must be described in terms of whole number 
values in the coordinate system. 

[In an attempt to linit problem complexity, structures such as turning vanes, 
noise attenuators, the test stand, and various other obstructions within the test cell 
were omitted from the problem. Structures such as the ejector tube and Coanda surface 
were greatly simplified, as were geometric shapes such as circles and ellipsoids. Where 
possible these more complicated shapes were replaced with flat, square, or rectangular 
shapes. 

The driving factor behind the grid dimension was storage space available in the 
PHOENICS code itself. Since all surfaces and known conditions must be identified 
with either a “PATCH” or “CONPOR” statement, the complex geometry of the 
Coanda cell required a quite lengthy “QI” file. The PITOENICS code, as currently 
installed at the Naval Postgraduate School, has a limit of 100 PATCH statements 
available. This dimension of 100 is divided among all PATCH and CONPOR 
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Figure 3.1 Gas Turbine Test Cell (Side View) 
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statements, with a CONPOR roughly equal to six PATCH statements in storage 
requirements. A grid definition finer than the one used in this problem or the inclusion 
of some of the omitted structures listed above, would exceed the storage dimension and 
prevent the successful use of the code. 

With these factors in mind, the problem was established using the following 


boundary conditions: 


® ambient pressure 101325 Mb 

* ambient temperature 295 K 

¢ mass flow (engine) 175 kg/sec 

® gas constant 286 m*/(sec*K) 
*c, 1004 m*/(sec“K) 


e jet blast temperature 2000 K 

eno slip” boundary condition on all surfaces 
These conditions were used in various places in the “QI” file to define the physical 
parameters of the cell during operation. The final form of the QI file used in this 
problem is presented in Appendix A. Specific statements in the QI! file are described in 
detail in the PHOENICS user manual [Ref. 4]. 

Initially the problem was solved without the use of the PHOENICS “turbulence 
niodel” in an attempt to achieve more rapid convergence to a solution. This laminar 
problem solved for the dependent variables pressure, enthalpy, and the three velocity 
components. Only after convergence was achieved was the more complex turbulence 
model invoked and the process repeated. The turbulent model solved for the five 
variables previously mentioned, plus kinetic energy and epsilon. Due to the extreme 
complexity of the model when the turbulence model is used, it was found that 
convergence was highly unlikely unless this laminar-turbulent routine was followed. 

The turbulence niodel used was the "KEMODL,” one of several models available 
through PHOENICS. By invoking KEMODL, the k-€ turbulence model is added to the 
problem. Two additional dependent variables, .urbulence kinetic energy (k) and 
turbulence kinetic energy dissipation rate (Cc), are solved. The turbulence calculations 
are governed by the following equations for turbulent viscosity and mixing length, 


respectively: 


= 
Hl 


1/2 
Cx 


— 


3/2), 
m Cok! 'e 
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where 


as = turbulent viscosity 

Cy = viscosity constant 

k = turbulence kinetic energy 

ra = turbulence mixing length 

eT) = dissipation constant 

€ = kinetic energy dissipation rate 


These equations, which have been developed to model a variety of turbulent flow 
situations, are thoroughly described in [Ref. 5: p.91]. 

The problem was solved using the MVS batch system of the IBM 3033 
mainframe at NPS. Due to the size of the Coanda problem, all runs of over 75 
consecutive sweeps were performed overnight using the appropriate class code. Due to 
the fact that PHOENICS writes over data stored in the DFO9 file during each run, 
DFO9 data was transferred to long term disk storage backup after each successful run 


in order to protect the vital restart data. 


C. MAJOR ASSUMPTIONS IN THE MODEL 

Several major assumptions were made in the modeling of this problem. The first 
major assumption was in the area of geometrical approximation. The two most 
noteworthy of these approximations were in the modeling of the Coanda surface and 
the ejector tube. The curved Coanda surface was modeled as a right angle, the sides of 
which connect the exit of the ejector tube and the bulkhead between the exhaust and 
the after secondary inlet. The ejector tube, which in reality is a duct of complex 
geometry with cooling inlets along its length, was replaced by a simpler rectangular 
cross-section duct with a porous surface. 

It was assumed that the deletion of the acoustic attenuators, turning vanes, and 
interior obstructions would have minimal effects on the general quality of the results. 
Also neglected were structures such as the supports for the Coanda surface and the test 
stand. These items, while altering the aerothermal conditions to a certain extent, are 
thought to play a minimal role in the overall performance of the test facility. 

The other assumption which merits discussion here concerns the question of grid 
cell size. As previously discussed, problem definition was driven by test cell geometry 


and limitations imposed by the dimensions within the code. The initial strategy of the 


ZI 


research was to solve the problem using a coarse grid and then enhance the solution 
through the use of a fine grid, time permitting. The limitations on problem size in the 
current PHHOENICS installation and complexity of the problem along with the limited 
time available for thesis research prevented the establishment of the fine grid model. 
While the grid cell size finally used in the problem is somewhat coarse, the assumption 
was made that an adequate solution could be obtained and later refined with a fine grid 
after the PATCH limitations of the PHOENICS code are expanded. 


LY. DISCUSSION OF THE RESULTS 


A. OBTAINING A FINAL SOLUTION 

As briefly mentioned in Chapter III, the solution of the problem was obtained by 
first using a laminar model, followed by a restart using the turbulent model. The 
laminar model, which solved for pressure, enthalpy, and the three velocity components, 
was run untl convergence was achieved. RELAX commands were used on pressure 
and the three velocities in an effort to force quicker convergence. When appropriate 
relaxation values were given, the problem converged in approximately 1000 sweeps. 

The turbulent model was then invoked bv 'nserting several additional commands 
into the QI] file. The significance of these commands is discussed in the PHOENICS 
literature. Lsing the five parameters (P,H,L,V,W) solved for in the laminar model as 
the data for a restart, the turbulent model was run. 

Initially, the increased complexity of the turbulent problem led to the divergence 
of the solution. The RELAX values had to be adjusted in order to insure convergence. 
The most difficult properties to contain in the turbulent model were enthalpy and 
kinetic energy. Even when the divergence was stopped, the problem was very slow to 
converge to a solution. The addition of two additional properties by the turbulent 
model, kinetic energy and epsilon, caused the computer time to increase drastically. 
Frequently, only 50 sweeps could be obtained in a 15 minute (NPS Class=G) time 
frame. Overnight runs of 60 minutes (NPS Class=J) would allow only 250 sweeps. 

The turbulent model converged in approximately 2500 sweeps after the restart 
from the laminar data. It should be noted however, that the sweeps were not taken in 
consecutive runs. Alterations to the QI file were required between runs in order to 
insure that PHOENICS continued to converge to a solution. Convergence was judged 
to have occurred when the residuals varied only slightly over a series of sweeps and the 
mass balance was satisfied. The final 50 sweeps taken in the turbulent model are 
presented in Appendix B. The output begins with information as to the grid spacing, 
values to be solved for, and properties defined in the QI file. The residuals for each 
sweep of the problem are then given, indicating the condition of the soluuon. The 
values within the test cell are given by “ZSLAB,” with each ZSLAB corresponding to a 


Z coordinate in the model. Values for the Y positions in each slab are labeled as 
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“TY =” while the X position is indicated by the column number. The centerline of the 
test cell is column | (IX=1). The tabular data should be used in association with 


Figures 3.1 and 3.2, which indicate the location of cach of the grid cells. 


B. COMPARISON OF THE SOLUTION WITH EXPECTED RESULTS 

In general, the solution of the problem was found to be in agreement with 
expected results. The direction of flow within the test cell as well as the temperature 
and pressure distributions were generally found to be within conceivable limits. 
Comparison of the PHOENICS output in Appendix B with Figures 2.2, 3.1 & 3.2 
shows the correlation between the model solution and the test cell concept. Figure 4.1 
shows values obtained at selected points and how they compare with the test cell 
concept presented in Figure 2.2. This comparison is for flow direction only, since no 


numerical field data was available for a quantitative con. aiison. 


REGION PRIMARY DIRECTION COMPARISON 
Primary Inlet W = -6.5 m/s agrees 
Forward Secondary WwW = 15.1 m/s agrees 
Iniect 

After Secondary V = 9.5 m/s disagrees 
Inlet 

Exhaust V= 28 mis agrees 
Test Room W= 14 m/s agrees 
Eyector Ibe W = 106 m/s agrees 


Figure 4.1 Comparison of Flow Direction between Results and Concept. 


As indicated by the data, the flow directions within the cell roughly coincide with 
the Coanda concept. The major difference lies in the after secondary inlet. The concept 
of the Coanda test cell is that cool secondary air in inducted through the after 
secondary inlet and mixes with the hot exhaust gases below the Coanda surface. This 
muxture then passes up the exhaust stack and out of the test cell. The results of the 
model show, however, that the flow in the after secondary inlet is out of the cell. 


Throughout the modeling of the test facility, this failure of the after secondary inlet has 


persisted. It is unclear at this time whether this anomaly is a problem with the test cell 
design or is a problem with the modeling of the facility. A closer investigation is 
certainly warranted in this area. 

While the actual magnitudes of the values within the cell are unknown, it is a 
safe assumption that the model failed in the region of the engine exhaust. The values of 
temperature and velocity (938 °K and 287 m:s respectively) are lower than those 
expected. These properties, which will be discussed in the next section, were found to 
have decreased significantly when the turbulent model was invoked. It should be noted 
however, that even though the values are low, the general flow pattern generated bv 
the solution is considered to be relatively accurate. 

The magnitudes of the properties in other regions of the test cell will undoubtedly 
change when the engine region problem is corrected. It was found that the solution of 
the entire problem changes drastically as the result of the altering of these key values. 
An attempt was made to dictate the precise conditions at the exhaust of the engine 
after the problem had converged, and the result was disastrous. Basically, the entire 
problem would have had to been completely reworked, beginning with the laminar 


solution... a prospect that a lack of adequate time prohibited. 


C. COMPARISON OF LAMINAR AND TURBULENT SOLUTIONS 

The successful completion of a laminar model was a necessary step in the 
solution of the problem. While the the problem is most certainly a turbulent one, the 
simpler laminar case was used to provide “ball-park” values for the restart of the 
turbulent model. As shown in Figure 4.2, the results of the laminar and turbulent 
models varied considerably. In both cases, the flow behavior generally followed the 
expected result, however the magnitudes of the properties were different. In the laminar 
model, the values around the engine were particularly good, with velocities and 
temperatures close to those expected from the F404 engine. Other values, however, did 
not seem to agree with the information gathered during the on-site inspection. The test 
cell rear wall temperature, for example, was found to be 1836 °k in the laminar model 
while test personnel reported the wall only warm to the touch immediately after engine 
shut-down. The turbulent model, on the other hand, gave believable values for the 
rear wall temperature while producing poor values in the vicinity of the engine exhaust. 

While the cause for these differences is not precisely known, it 1s thought to arise 
from the large grid cell sizes taken in the model. Due to the nature of turbulence, 


particularly in the extremely violent conditions at the exhaust of the engine and ejector 


LOCATION PROPERTY LAMINAR TURBULENT 
engine inlet T (°K) 3125 2o 
ae -1202 MPa -1453 MPa 
velocity( W) 175 m/s 175 m/s 
engine exhaust T (°K) 2000 K 938 K 


AP 4377 MPa 12100 MPa 
velocity( W) 745 m/s 207 TS 
test cell rear wall T (°K) 1836 K 482 K 
test cell exhaust T (°K) 361 K 447 K 
velocity( V) 109 m/s 20 mis 





Figure 4.2 Values at Key Locations in the Test Cell. 


tube, the scale of the grid must be made smaller in order to capture the phenomenon. 
In this model, the dimension of the first grid cell after the engine exhaust was I meter 
in the Z direction. A more accurate solution might have arisen had this single grid cell 
been divided into several smaller cells. This concept of smaller grid sizes would 
probably yield better results at not only the engine region, but at any region of extreme 
turbulence, pressure, or temperature. While PHOENICS gives general guidance as to 
the final grid size, the establishment of the final grid is basically a “trial and error” 
procedure in an attempt to find a grid independent solution [Ref. 4: p.8.11]. 

The solution presented is considered to be a good starting point for further 
investigation of the Coanda test cell. While the quality of some of the magnitudes in 
the solution are questionable, the basic concept of the test cell seems to be well 
represented and therefore indicates the ability of PHOENICS to solve the problem. 
The obviously erroneous data which was generated in the engine exhaust region is 
considered to be a local problem within the cell solution, and is not indicative of the 
failure of the code to obtain results elsewhere in the problem. While the poor results in 
the engine region do affect the results “downstream”, the properties “upstream” seem to 
have remained virtually unaffected. Even though some of the data in the laminar 
solution initially appears more accurate, the true values of the solution will 
undoubtedly be derived from the turbulent solution after the correction of the engine 


exhaust problem. 
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V. SUMMARY AND RECOMMENDATIONS 


A. ASUMMARY OF THE PROJECT 

In an effort to analyze the performance and determine the design effectiveness of 
Gas Turbine Test Cell #1, Lemoore NAS, California, a numerical analysis of the 
Operating parameters was performed using the Parabolic, Hyperbolic, or Elliptic 
Numerical Integration Code Series (PHOENICS). The project involved the on-site 
inspection of the facility, modeling of the test cell, and the utilization of the code to 
determine the physical conditions within the cell during the full power run-up of a GE 
F404 gas turbine engine. The analysis was first performed using a laminar model of the 
test cell, followed by the “restart” of the code with a turbulence model invoked. Data 
was generated in tabular format. This data, in conjunction with the grid model 
diagram, gives the user insight into the aerothermal characteristics of the Coanda Test 
Cell 

The results gained through PHOENICS were generally good, with the model 
roughly replicating the expected results. Obvious errors did occur, however, in the 
region of the engine exhaust and ejector tube following the introduction of the 
turbulence model. This problem is thought to have resulted from a coarse grid cell size 
in regions of extreme turbulence, temperature and/or pressure. 

The flow within the model differed from the expected result in the after secondary 
inlet. Instead of indicating the induction of secondary cooling air through the inlet, the 
model indicated that air was flowing the opposite direction and exiting the cell. 
Currently, it is not known if this difference is due to modeling error or a design error in 


the test cell. This matter warrants further investigation. 


B. RECOMMENDATIONS FOR FURTHER RESEARCH 

From the groundwork laid by this thesis, it is apparent that the PHOENICS code 
is a useful tool in the solution of the problem. Several changes to both the model and 
the installation of PHOENICS at NPS should be made prior to continuing work on the 
numerical simulation of Lemoore Test Cell #1. It 1s believed that the model presented 


here will produce better results if the following alterations are made. 


aah 


1. Improvements to PHOENICS at NPS 

Overall, PHOENICS performed well during the course of the research. It 
should be noted however, that some limitations in the code as it is currently installed at 
NPS did cause some difficulties. The dimensions limiting the maximum allowable 
number of PATCH and CONPOR statements should be increased from the current 
value of 100 to at least 250. This change would allow a more accurate definition of 
complex problems such as the Coanda test cell. With the option of using more PATCH 
and CONPOR statements, many more boundary conditions could be established within 
a relatively coarse grid or a much finer grid with associated boundary conditions could 
be used. Either option would undoubtedly lead to an improved solution. 

In addition to the dimension change, the PEHIOTON graphics option of 
PHOENICS should be made available. The ability to produce a graphic representation 
of the physical conditions within the problem would be an invaluable asset. 

2. Improvements to the Coanda Model 

The most obvious improvements to the Coanda model are in the area of grid 
size. Grid cells should be made much smaller in the engine exhaust region, the ejector 
tube, all inlets and around the Coanda surface. This reduction could be achieved in 
several ways, the most apparent of which is the redefinition of the problem with a 
different grid spacing. The difficultv in problem redefinition lies in the fact that the 
establishment of a revised grid system is very complicated and time consuming. 
Another approach, which might allow for the current grid to be used, would be 
through the use of “sub-models” over regions of special interest. The results of these 
sub-models could then be reinserted into the original test cell model as boundary 
conditions. In addition to eliminating the need for the complete rework of the test cell 
problem, this method may reduce computer run times and bypass the need for 
PHOENICS code dimension alteration. 

Additional field data would also greatly enhance the ability of the model to 
more accurately determine the conditions within the test cell. Values such as the mass 
flow at each inlet and the exhaust, pressures and temperatures at critical points within 
the cell, more specifics about the engine performance, and information on areas of 
recirculation could be inserted as boundary conditions into the current model. The 
establishment of these additional boundary conditions would help to insure that 


PHOENITCS is converging to the most realistic solution possible. 
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©.) CONCLUSION 

This thesis produced results which represent a foundation on which the 
completion of the project can be based. The results obtained indicate the successful 
modeling of the gas turbine test facility and suggest ways in which PHOENICS may be 
more effectively used in the completion of this problem and in the solution of others. 
With relatively minor changes to the current model of the test facility, accurate values 
for the aerothermal conditions within the cell may be obtained. 

Judging form the results obtained to date, the Coanda test cell performs as it was 
intended to perform. Questions arise however, concerning the effectiveness of the 
Coanda surface and the after secondary inlet. Closer examination of these regions 


should reveal the actual conditions present during test cell operation. 


Pes, 


APPENDIX A 
THE PHOENICS Q!1 FILE 


TALK=F ;RUN(1,1) 


KEKKKKKKRKKKKKKKKKKRKRKKRRKRKKRAKKRKRRKRRRKKRRKRRRKRRKRRKRRKRKRKRRRRKKKRKKRKRARRARKAREK 
KREAKRKRKRKRRARKRKAKRKRKRRRKRRRRKRKRKKKRKRKRKRKRKRKRRRRRERRRERRRRRRRRRRRRRRRKRRRRRRRRKI 


AK* ARK 
RRR LEMOORE TEST CELL #1 (COANDA TYPE) RRR 
AKK AKK 
RRR COARSE GRID AK 
RR TURBULENT MODEL RAK 
AKKX RR 


RAARKEKKKRRKKRKKRKKKRKRKRKEKRKRRKRRRRRRRRKRRERKRKRRRRKRRRRRKRRRRRRRRRRRRRRRKRRRARERRKI 
KREAKKKRAKRKKKKKRKRKARKRKRKKRKRKKRKRRKRKRKRRKRRRRRRRRRRKRRRRKRRRKRRRRKRRRRKRERKERRKK 
* 


GROUPL: TITLE 
TEXT(COANDA RUN # 3 ; COARSE GRID) 
REAL (WJET MJET TJET, PJET ,RHOJET,HJET .GASCON, CSUBP, TAMB, WIN) 


REAL ( KEINJ,EPINJ,KEINA,EPINA 
GROUP2: STEADY/TRANSIENT 


STEADY =T 

GROUP3: X DIRECTION GRID 
NX= 6 
XFRAC(1)=-1. ; XFRAC(2)=0.5 
XFRAC(3)= 1. ; XFRAC(4)=1-41 
XFRAC(5)= 2. ; XFRAC(6)=0.685 
XFRAC(7)= 1. ; XFRAC(8)=0.305 
XFRAC(9)= 1. ; XFRAC(10)=1.0 

GROUP4: Y DIRECTION GRID 
NY = 16 
YFRAC(1)=-1. ; YFRAC(2)= 0.15 
YFRAC(3)= 1. ; YFRAC(4)= 0.815 
YFRAC(5)= 1. ; YFRAC(6)= 1.00 
YFRAC(7)= 1. ; YFRAC(8)= 0.820 
YFRAC(9)= 1. ; YFRAC(10) = 1.19 
YFRAC(11)=1. ; YFRAC(12) = .305 
YFRAC(13)=1. ; YFRAC(14) = 1.52 
YFRAC(15)=1. ; YFRAC(16) = 0.305 
YFRAC(17)=1. ; YFRAC(18) = 0.305 
YFRAC(19)=1. ; YFRAC(20) = 0.835 
YFRAC(21)=1. ; YFRAC(22) = 0.460 
YFRAC(23)=1. ; YFRAC(24) = .305 
YFRAC(25)=2. ; YFRAC(26) = 2.78 
YFRAC(27)=1. ; YFRAC(28) = 3.66 

; YFRAC(30) = 1.00 


GROUPS: 2 DIRECTION GRID 


NZ = 19 

ZERAC(1) = -2. ; ZERAG\ 2) —sivez 

Z2FRAC(3) = 2. ; ZFRAC(4) = 2.44 

ZERAC(5) = 2. 7; ZERRG(Gpe sls 

ZFRAC(7) = 1. ; ZFRAC(8) = 0.305 
ZERAC(9) «=m 1. ; CRRAC( 10) — sire 

ZERAC(II)= 1. ; 2ERRCtIZ)— 07205 
CFRAC(I3)= 1. ; ZPRAC(WI4)— 33505 
ZFRAC{IS)= 1. ; ZFRA@OTSG)= 1.00 
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ZFRAC(17)= 1. ; ZFRAC(18)= 0.305 
ZEREC(19)= 1. ; ZPRAC(20)= 1.22 
@armerct)= 1, ; ZFRAC(22)= 0.305 
ZFRAC(23)= 1. ; ZFRAC(24)= 3.68 
@pemeves/= 1. ; ZFRAC(26)= 0.457 
@eene(zZ/)= 1. ; ZERAC(28)= 3.0 
ZFRAC(Z9)= 1. ; ZFRAC(30)= 0.457 
2ERAC(31)= 1. ; ZFRAC(32)= 3.0 


GROUP6: BODY FITTED COORDINATES 
GROUP? : VARIABLE SOLUTION, STORAGE & NAME 


Somve Pl ,H1,U1,V1,W1) 
SeeuIN(P1,Y,Y,Y,N,N,N 
Seen ,Y,Y 
SoEUrN(UL,Y,Y 
SsemurTi vl, Y,Y 
SeLUri(wl,Y,Y 
STORE (RHO1, TM 


P 
GROUP8: T 
GROUPS: P 
295. 
2000. 


N 
,N 
/N 
ei 


AAAs 


c 
‘ 
, 
c 


sg 
N 
N 
N 
1 
E 


RMS IN EQUATIONS 
HYSICAL PROPERTIES OF MEDIUM 


TAMB 
TIT 
Bjer = 0. 
Bier = 175. 
PRESSO = 101325. 
feos = 0. 
GASCON = 286. 
RHO1B = 1./GASCON 
RHOJET = (PRESSO + PJET)* RHO1B/TJET 
WJET = MJET/RHOJET 
THPLA = TINY 
CSUBP = 1004. 
TMP1B = 1./CSUBP 
HJET = CSUBP * TJET 
TMP1 = GRND2 
RHO1 = GRNDS 
DRH1DP = GRNDS 
ENUL = 1.E-05 
ENUT = 200. * ENUL 
WIN= MJET/1.0 


TURMOD (KEMODL) 

STORE (VIST) 

KEINA = 0.5 * (0.005*WJET)**2 

EPINA = 0.09 * KEINA **1.5/(0.1) 
Rn 


Grouelg: INTBBPRASE PROCESSES 
GROUPILI: INITIAL VALUES 


eee) = 1.0 
BPLLNLT EE KEINA 
FIINIT(EP EPINA 
xn 


Poet INIVAL,1,2,2,4,8,19,1,1) 
Peer wees ,0.0,0.2*WJET 
Beetowihiit,m.,0.0,CSUBP*0.2*IJET) 
eet T ee no Re ot 

Pee (Vie BP, 0.0,EPINJI*0.2 


**X INSURE THAT RESTART VARIABLES MATCH THOSE CURRENTLY STORED *** 
an Pitre PrO9 FILE... 


RESTRT (ALL) 
SECU +2: ADJUSTMENT TO FLUXES 
GROUPL3: BOUNDARY CONDITIONS 


* 
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KRARKKAAKKKKKKKAKAKAKKAKAKKAKKKKKAKKKKKKAKKAKKAKKKKKKAKKAKKKKKAKAKKKKAKKKAKRK 


SaRKK KRAKK 
RARKK ESTABLISHMENT OF PHYSICAL DIMENSIONS RAKKK 
KAR RARKK 


RAKRKKKKAKKAKKKKKKKKKKKKKAKKAKKKKKRKKAKKRKRKKRKKRKKAKRKKKAKAKKKRKKKARRKKRKKAKRKKRKRKAKRK 
x 


5 kkk TEST CELL DECK *** 


PATCHCDECK ,SWALL, 15415, [elo ole) 
COVAL(TDECK , Ul, GRND2Z, 0.0 

COVAL (DECK ,W1,GRND2,0.0 

COVAL BE Can en 

COVAL DECK, ,EP,GRND2,GRND2 


x xX TEST CELL FRONT WALL *** 


PATCH(FRONT, LWALL,1,4,1,11 eee) 
COVAL (FRONT, U1,GRND2Z,0.0 

COVAL( FRONT ,V1,GRND2,0.0 

COVAL Be ee eae el 

COVAL FRONT ,EP,GRND2Z,GRNDZ 


* *xk*k PRIMARY INLET TOP *** 


PATCH(PRITOP ,NWALL,1,4,11,11,1,4,1,1) 
COVAL(PRITOP,U1,GRND2,0.0 

COVAL(PRITOP ,W1,GRND2,0.0 

PATCH( PRIROOF ,SWALL,1,4,13,13,1,4,1,1) 
COVAL(PRIROOF,U1,GRND2,0.0 
COVAL(PRIROOF ,W1,GRND2,0.0 
PATCH(PRIEND, LWALL,1,4,12 
COVAL PRS Dut Giere a 
COVAL(PRIEND,V1,GRND2,0.0 


x *xkX PRIMARY INLET DIVIDER *** 


PATCHCPIDIOP, SWALE, bea oaee2 4,1,1) 
COVALPIDIOP, UL 7GRNBZ, 020 
COVAECPIDIOP Wl GaNbZ Ono 
PATCHVELDBOT ,NWALE I 4, (ee. 4, 1,1) 
COVAL( PIDBOT , U1, GRNib2, 0.0 
COVALCEIDBOT Wl Ganz, 0.0 
PATCH( PIDLFT ,HWALL,1,4,8, 
COVAECPIDEET , UL TGRiDZ . 0.0 
GCOVALE(PIDLET ,VI{GRNDZ, 0.0 
PATCH( PIDRGT,LWALL,1,4,8 
4 


eZ), 5 2) 


COVAL( PIDRGT ,V1,GRND2 
PATCH( PRIINB,SWALL,1, 
COVAL( PRIINB,U1,GRND2 
COVAL PRIINB,W1,GRND2 


. *x*xX TEST ROOM OVERHEAD *** 


PATCHUDROVHD ,NWAbIn I 4,5 ,5) Se 11) 
COVAL(TROVHD ,U1,GRND2,0.0 
COVAL(TROVHD ,W1,GRND2,0.0 

PATCH( TRLFT,HWALL,1,4,6,6,2,2,1,1) 
COVAL(TRLFT,U1,GRND2,0.0 

COVAL TRLFT,V1,GRND2,0.0 


a *kX PRIMARY/FORWARD SECONDARY INLET DIVIDER *** 
PATCH INBEOGM1 , LWAUEy ade o> i Se Ss ee 
COVAL( INFLOW] ,W1,FIXVAL,-6.5 
Pane ENE LOW2 , LWAUEpel 47 7 oo 
COVAL INFLOW2 ,W1,FIXVAL,-6.5) 


BAT CMBRSIDL - HWARIaL 2,7 12 ,6)6)1,1) 
COVAL(PFSIDL,U1,GRND2,0.0) 
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COVAL 
PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 


PaLcH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 
Pee CH 
COVAL 
COVAL 


PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 

x 


x 


* 


PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 
PATCH 
COVAL 


COVAL(TRURER, vi GRND2, 


x 
PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 
PATCH 
COVAL 


COVAL(TRLRER, V1,GRND2,0. 


E@erDE, V1l,GRND2Z,0.0) 
RostDt ovack, 1, 4, Per 37, 7,1,1) 
Peso. ,UL, GRND2, 0.0 
eowDr, Wi, oe O26 
PeSiDR, LWALL, ime le, o,6,L,1) 
weoeDR, Ul, GRND2, ce 0 
PeolpR, VL;7GRNDZ,0.0 


*x** FORWARD SECONDARY INLET TOP *** 


FSLLIP,HWALL,1,4,9 
Hope P UL, GRND2Z,O 
HoeLtP V1, GRNDZ,O. 
FoSO! , SWALL, 14,7, 
Foo. Ul, GRND2, 0.0 
FopOr WL,GRNDZ,O.0 
FSOVHD, NWALL, 1,4 : 
PSOVHD , Ud ,GRND2, 
FSOvED, Wl, GRNDZ, 
POLLY, SWALL, M 4,1 
BothiP Ui ,GRND2,0. 
Motu le Wi eanib2, 0. 
PSRELP , LWALL, Pea 
FORLIP , Ul, GRND2, e. 
HSREEP,Vi,GRND2Z,0. 
FSMOF ; SWALL , t, 2,20 
FSTOP,U1,GRND2,0. 0 
FSTOP ,W1,GRNDZ,0.0 


OMS Re sl pl) 


} 


0 
0 
a 


a) 


oo 


io kL) 
Oro, LO, 1, L> 


OrOe 4), 1,1) 


*x*X FORWARD SECONDARY TURNING BULKHEAD *** 


Tebeen,AWALL,1,4,6,8,13,13,1,1) 
TERESR , Ul, GRNDZ ,0.0 
TBREAR,V1,GRNDZ, 0 

TEOUT, LWALL , eo, 15, 15,11) 
EBOUT Ul, GRND2, 6! 0 

TBOUT ,V1,GRND2,0.0 

TBRRR, LWALE,1,4,5,6,13,13,1,1) 
PeeRNT ,UL,GRNDZ,O.0 


LTBFRNT ,V1,GRNDZ,0.0 


Riot ROOM ART BULKHEAD *~** 
UPPER 


TRUPFT ,HWALL,1,4,4 
TRURET U1, GRND2, o- 

TRUPFT ,V1,GRND2, Oe 
TRUBOT ,NWALL,1, 4 : 
EReUBOL, Ul, GRND2, 
BeupOr, Wl, ,GRND2, 
TRURER, LWALL, i, 4,4 
Proueen, U1, GRND2, 


OO- 


LOWER 


TRE | HWALL 14,1 
TRLFT, Ul ,GRND2,0. 
Meert Vi, GRNDZ ,O. 
ER GLP ; SWALL, 1,4,3 
OSErP ,U1 ,GRND2,0. 
Droeer  Vl1,GRND2 ,O. 
TRLRER, LWALL, 1,4, 
TRoeeR , Ul ,GRND2,0. 


leaky tpl) 


| 
mere 2 PZ ll) 
$ 
0 
ee, 1313515) 
0 
0 


33 


PAT GE 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 
PAECH 
COVAL 


Srpe 


VIRSET ,HWALL {274s Lee 
Tent, ol ,GRND2,0. Q 
TROPT , Vl ;GRNDZ O28 
Role, EWALL, 1 17353 pees) 


,o 
RSE Vt ,GRND2,0. a 
TRSTP,W1,GRNDZ 0.0 
TRSRER, LWALL, Zs of 

0 


15,3, 13 eee) 
TRORER,UL /GRND2 a 


COVAL( TRSRER, V1 "GRND2, 


xk ENGINE **% 


* 
CONPOR(0.0,CELL,1,1,3,3,9,10) 


PATCH 
COVAL 
COVAL 
COVAL 


oan 
COVAL 
COVAL 
COVAL 
COVAL 
COVAL 

x 


* 


x 


PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 


x 


SEEN FL; FIXFLU, -MJE 
JETIN,H1 ,ONLYMS , SAME ) 
JETIN,W1,FIXVAL WIN) 


JETOUT ,LOW 171355 eel 
JETOUT , Pl, FLXAREU ioe 

JETOUT ,W1,ONLYMS ,WJET 
JETOUL, Hl _OCULYUS , Muir 

fete eee ey 
JETOUT EP sOneyis  Sring : 


PStiN HIGH tl 3,35.6, nie yl) 


xk EJECTOR ROOM OVERHEAD *** 


EROVHD UL, GRAD2,6.0)" 0 
EROVHD ,U1,GRND2,0. 

EROVHD,W1, ,GRND2, Ole 

ERROOF, SWALL , i, 4, Te "4, oo lord; 1) 
ERROOE, Ul ,GRND2, 0.0 

ERROOF ,W1,GRND2,0.0 


kkk EJECTOR TUBE *** 


CONPOR(O-.S RASTA, is,57027 lo) 


GCNEOR( 01S NORTH, 1403, 3, 


Ze aes 


CONPOR(0.5,SOUTH,1,1,3,3,12,16 
n 


* 


* 


PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 


* 


PATCH 
COVAL 
COVAL 
PATCH 
COVAL 
COVAL 
PATCH 
COVAL 


COVAL(ERLRER, V1,GRND2 


xeX EJECTOR ROOM APT BULKHEAD *** 
UPPER 


ERUFT ,HWALL,1,4,4, 
BaUeT UL, GRND2, 0.0 
BPRURT Vil Gaubz 070 
ERUBOT, NWALL, 1,4,3,3 be, en 
BRUBOT ,Ul ,GRND2,0. 0 
ERUBOT,W1,GRND2,0. 0 


LOWER 


ERLFET,HWALL,1,4,1 
ERLFT, UL ,GRND2, 0. 
ERLFT,V1 ,GRND2, O. 
ERLTOP /SWALL, 1,4 

4 


pLoue onal) 


~~” U1 


, Lo phos.) 


, Loe lerl, 1) 
ERLIOP UL, GRND2 
ERLTOP, Wl GRiDZ 
ERLRER, LWALL, 3 

ERLRER,U1 ,GRND2 


oO O-~ 


SIDE 
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PomenvsksrL,HWALL,2,4,3 
GeVenyenor.l,Ul, GRND2, 


i pees 2) 
COVAL(ERSFT,V1,GRND2,0. "04 
Is oe 


Peter ERSLIP, EWALL, i. poplar l6>.1,1) 
GOvVALCERSLIP,V1, GRND2, 518) 
COVA ERSLIP (WL, aan a) 
PATCH ERSRER, LWALL, o 4°3°3,17,17,1,1) 


COVAL ERSRER,U1,GRND2, 0.0 
COVAL ERSRER,V1,GRND2,0.0 
x 


a *x*xX AFTER SECONDARY FORWARD BULKHEAD *** 


PATCH(ASINR,LWALL,1,4, 5, aie 1 
COVAL(ASINR,U1,GRND2,0. 
COVAL(ASINR,V1,GRND2.0. 5 
PATCH(ASOUTR,HWALL,1,4,7,14,15,15,1,1) 
COVAL(ASOUTR,U1,GRND2,0. 8 
COVAL(ASOUTR,V1,GRND2.0.0 
PATCH(ASTOP,SWALL,1, 4, 15,15,16,16,1,1) 
COVAL(ASTOP,U1,GRND2, a 

COVAL ASTOP WI GRND2.0 

x 


| xk AFTER SECONDARY/EXHAUST DIVIDER *** 
x 


PATCH(ASEDFT,HWALL,1,4,10,15,17,17,1,1) 
COVAL(ASEDFT,U1,GRND2,0.0 
COVAL(ASEDFT,V1,GRND2,0.0 
PATCH(ASEDTP ,SWALL,1,4,16,16,18,18,1,1) 
COVAL(ASEDTP ,U1,GRND2,0.0 
COVAL(ASEDTP ,W1,GRND2,0.0 
PATCH(ASEDRR ,LWALL,1,4,10,15,19,19,1,1) 
COVAL BSEDRR U1, GRNDZ, 0.0} 
COVAL(ASEDRR,V1,GRND2,0.0 
COVAL(ASEDRR,KE,GRND2,GRND2 
COVAL(ASEDRR,EP,GRND2,GRND2 
PATCH(ASEDBT ,NWALL,1,4,9,9,18,18,1,1) 
COVAL(ASEDBT ,U1,GRND2,0.0 
COVAL(ASEDBT ,W1,GRND2,0.0 

x 


p xk TEST CELL REAR WALL *** 


PATCH(REARWL ,HWALL,1,4,1,15,19,19,1,1) 
COVAL(REARWL,U1,GRND2,0.0 
COVAL(REARWL,V1,GRND2,0.0 

COVAL pte ge GONDZ . GRNDZ) 
COVAL(REARWL ,EP,GRND2,GRND2 


Be *x*xX TEST CELL INSIDE WALL *** 


Bement o lpr EWALL,4,4,1,11,1,7,1,1) 
COVAL(SIDE1,V1,GRND2,0. 
COVAL(SIDE1, Wl, GRND2,0. 
Peach SlDEZ, EWALL, 4, ete) Brot) 
COVAL SIDE2,V1,GRND2,0. 0 
COVAL(SIDE2 ,W1,GRND2Z,0.0 


OO 
ee 


Papers sIfh2, EWALL,4,4,1,8,10,14,1,1) 
COVAL(SIDE3,V1,GRND2,0.0 
COVAL(SIDE3,W1,GRND2,0.0 
PATCH(SIDE4,EWALL,4,4,1,5,15,15,1,1) 


COVAL(SIDE4,V1,GRND2, 
COVAL(SIDE4 ,W1,GRND2, 
Pewentvolves,EWALL,4, 4,1 ; 
COVAL SIDE5,V1,GRND2,0. i 


eed 
- OO 


Ae16,17,1,1) 


$) 
1 
COVAL(SIDE5,W1,GRND2,0. 
PATCH({SIDE6,EWALL,4,4,1,15,18,19,1,1) 
COVAL(SIDE6 ,V1,GRND2,0. ah 
COVAL(SIDEG, Wile Gripe 0.0 


: 
*xzX TEST CELL OUTSIDE WALL *** 
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PATCH (OU WWALL) Gpapie, 2 | ly elie) 
COVAL(OUT1,V1,GRND2,0.0 
COVAL(OUTI1, Wl, GRND2,0.0 
PATCH(OUTZ, /WWALL , BF 6, lO. 879. 11) 
COVAL OUT2,V1,GRND2, 0.0 
COVAL OUT2,W1,GRND2,0.0 


PATCH(OUT3 ,WWALL,6,6,1,9,10,14,1,1) 
COVAL Se ace aoa ay 
COVAL(OUT3 ,W1,GRND2Z,0.0 
PATCH(OUT4 ,WWALL,6,6,1,6,15,15,1,1) 


COVAL(OUT4,V1 /GRND2,0. 0 
COVAL(OUT4,W1,GRND2Z,0.0 
PATCH( OUTS /WHALL, 6,6,1, 14516 tiie) 
COVAL(OUTS, ce GRND2,0.0 
COVAL OUTS, W1,GRND2,0.0 
PATCH(OUT6, WWALL, 6 Oy, lo ale ele) 
COVAL(OUT6, V1 /GRND2,0. 0 
COVAL(OUT6,W1,GRNDZ,0.0 
PATCH WIOPL, SWALL , ST Ogres 
COVAL(WTOP1,UL1 ,GRND2,0.0 
COVAL(WTOP1,W1,GRND2Z,0.0 
PATCHUMIOR 2. SWALL, 575, 11, 869. 1ea9) 
COVAL(WTOP2Z,U1 ,GRND2,0. 0 
COVAL WTOP2,W1,GRND2 0 
PATCH(WTOP3 ,SWALL,5, 5,10 
COVAL WIOP3,U1,GRND2,0.0 
COVAL(WTOP3 ,W1,GRNDZ ,0.0 
PATCH(WTOP4,SWALL,5,5,7 
COVAL WIOP4 ,U1,GRND2,0. 
COVAL(WTOP4,W1,GRND2Z,0. 
PALCHMILOF 5, SWALE: 5 on 
5 


3) ale) 


—— 


LO; LO Lae) 
PLS, ope 


See el 
COVAL(WTOPS ,U1,GRND2 
COVAL(WTOPS ,W1,GRND2 
PATCH(WTOP6 ,SWALL,5,5, 
COVAL(WTOP6 ,U1,GRND2,0. 
COVAL (WTOP6 ,W1 ,GRND2 ,0. 


oy Sloe) 


. *k*X APPROXIMATION OF THE COANDA SURFACE *** 


CONPOR(0.0,NORTH,1,2,3,3,17,17) 
CONPOR(0.0,LOW,1,2,4,9,18,18) 
CONPOR(0.0,EAST,2,2,3,3,17, 12 
CONPOR 0.0,EAST,2,2,4,9,18,18 


. *XX BOUNDARIES TO AMBIENT *** 


* 


PAT CHU Ste NORTH 21, 4° Gy lon by Loe, 
COVALCSKY P10. i 

GOVAL (aig al, ONLYNS , CSUBP*TAMB ) 
COVAL(SKY,KE,ONLYMS, ern 

COVAL SKY, EP ,ONLYMS , EPINA 

x 


* 


PALGH (SOK ieee | 656,), 16,1. loei) 
COVAL(SKY1 Pl, Only 0.0 

COVAL(SKY1 /H1 /ONLYUS , CSUBP*TAMB ) 
COVAL(SKY1, PRE; ,ONLYMS ,KEINA 

COVAL SKY1,EP,ONLYHS , EPINA 


PATCHGSKVEEET, LOW.) G13 167i) 
COVE Siayweee l, Pi oO. 1 0.8 

COVAL(SWEEt nl Oilers CSUBP“TAME) 
COVAL Sedeer ce CN eeT) 

COVAL Swroee ll, EP OND Yio SPINA 


COUa eRe a OCOn tee aa 
COMBE(SKYRI, P1L,0.1,0:.0) 


COVAL(SKYRT,H1,ONLYMS , CSUBP*TAMB ) 
COVAL(SKYRT ,KE,ONLYMS , KEINA 
COVELL SKY RT, EP,ONLYMS, EPINA 
nw 


x 

GROUP14: DOWNSTREAM PRESSURE 

GROUP15: TERMINATION FOR INNER ITERATIONS 
x 


poree: = S01; LSWEEP = 550 
GROUP16: TERMINATION FOR SWEEPS & OUTER ITERATIONS 
GROUP17: UNDER RELAXATION 

Rm@bex(Pl,FALSDT,1.E-3 

Reb Ul, FALSDT,7.E=-2 

BELAX(V1,FALSDI,7.E-2 

BELAA(W1,FALSDT,7.E-2 


GROUP18: LIMITS ON VARIABLES 
wx 


VARMIN(P1)= -2.5E+05 
VARMAX(P1)= 2.5E+05 
VARMIN(V1)==-500 


VARMAA(V1)=500. 
VARMIN(W1)=-200. 
VARMAX(W1)=1200. 
VARMAZ = 
VARMIN(TMP1)=295. 
VARMAX(H1)= 2.008E+06 
VARMIN(H1)= 2.962E+05 


VARMIN(U1)=-200. 
VARMAX(U1)=200. 
VARMIN RHOL)=0. 
VARMAX(RHO1)=2. 
VARMAX (KE )=1.E+06 
VARMIN/ KE )=1 .E-08 
VARIIAX| BP =1 E+06 
VARMIN(EP)=1.E-08 
Kw 


GEOUrRILI: SeECIAL CALLS TO GROUND 

GROUFZO: PREGEMINARY PRINTOUT 
ECHO=T 

eeOQueZi: FIELDSPRINTOUT 

GROUP 22: MONITOR PRINTOUT 


TSTSWP=1 

IXMON = 2; IYMON = 5 ;: IZMON = 5 
GROUP23 : RBELD PRINTOUT & PLOTS 

ITABL = 3; NPLT = 15;IPROF =3 


GROUP24: PREPARATION FOR CONTINUATION RUNS 
OP 


oH 


APPENDIX B 
PHOENICS OUTPUT 


CCCC HHH PHOENTGS VERSION gis5) 03> obe? -i2co 
CCCCCCCC Hin (C) COPYRIGHT 1984 
CCCCCCC HAHREHAAAE CONCENTRATION HEAT AND MOMENTUM LTD 
CCCCCCC HHHBRRAeHnAn ALL RIGHTS ResenveED. 
CCCCCC HHHHHHHHHHHHH CHAM LTD, BAKERY HOUSE, 40 HIGH ST 
CCCCCCC HHHHHHHHHHHH WIMBLEDON, LONDON, SW19 SAU 
CCCCCCC BEAARnAinH TEL: 01-947-7/651 -.¢GmEEX: S2d517 
CECCCCCC annnn FACSIMILE: 0158/3 -e-e7 
CCC * nnn THE OPTION LEVEL Is od 


THIS CODE MAY ONLY BE USED UNDER THE TERMS AND CONDITIONS 
OF A LICENCE AGREEMENT WITH CHAM LTD. 


REPLICATION OF THIS CODE IS PROHIBITED UNLESS 
SPECIFICALLY AUTHORISED IN WRITING BY CHAM LTD. 


KRAKKKKKKKKKKKKKKKKKKKRKEKREKRKRKRKRRERKRKRKRRKRKRKRRRRKRKRKRKRKRKRKRRKRRKREKRKEKRI 


GREX1 OF 15/07/86 HAS BEEN CALLED 
FORMATTED SATLIT DATA READ FROM DF10 wh OR IRUN= 1 
pied eioticias STORAGE INFORMATION ------%* 
F "DIMNGN 200000 OCCUPIED= 51094 ESTIMATED MINIMUM DIMNSN= 15094 
DEP VRBL= 25536 OLD VRBL= 0 3D COEFF= 9024 3D DVDPS= 5472 


INITIAL FIELDS READ FROM DFO9 
KKAAKKKKKKKAKKKKKKKKKKKKKKAKEKKKAKKKKKKKEKKAKKKKKKERKKAKKRKKKK AK K 


GROUP 1. RUN TITLE & NUMBER 
KKKKKAKAAKKKKKKKAKKKEKKRKKKERRKEKKAKEKKKRRKAEKERKRKKKEREKAKKKERRKK 


RARKKKKKKKKRKKRKRKKKKRKKKKRKKRKKKKKKRKRKRKKRRKRKKRKRKRKRKRKRKRKRKKRERKRKKEKRKRRKKKKEKKE 


TEXT(COANDA RUN # 3 ; COARSE GRID ) 


KAEKKKKKRKKKKKRKRKKKRKEKRKKKRRKKKRRKRRKKRRRRKKRKRRRKKRKRRKEKRRRRKRRRRKEKRKKKKRKK 
RAEKEKKAKKRKRKKRKRKKKKKKKKRKRKRKRERRKRKRKRKKRERRKRKRKRRRRRRKRKKKRRERKRERKRKKKKKKKI 


IRUNN = i 
KAKKKKKKKKKKKKKKKKKKKKKKKKKKKRKEKEKKKKKKKKREKKAKKKKEKKKKKKKKKKK 
GROUP 2. TRANSIENCE 

STEADY c 

AKER KRAR RARA AERA RAR AAT IR IR Re EOE ee er 


GROUP 3. X-DIRECTION GRID SPACING 
T 


CARTES = 

a = 6 
AULAST = 1.000E+00 

METHOD OF PAIRS USED FOR GRID SETTING 
AFRAC 1) = -1.000E+00 ;XFRAC 2) = 5.000E-t2 
AFRAC 3) = 1.Q000E+OO ;XFRAC 4) = 1.410E+00 
AFRAC 5) = 2.000E+00 ;XFRAC 6) = 6.850E-01 
AFRAC 7) = 1.000E+00 ;XFRAC 8) = 3.050E-01 
AFRAC 2 OOE+00 :XFR 


120 AC 10 -OOOE+00 
RKAKKRKKKK KARR HHRA RAE KEI RERRER ARAMARK KARAT RRR IRA RRR 
GROUP 4. Y- Cee en GRID SPACING 


1.000E+00 
METHOD OF PAIRS USED FOR GRID SETTING 
YFRAC 1) = ~-1.,000E+00)- VERAC 2) = 1,500E=00 
YFRAC 3) = 1.000E+00 ;YFRAC 4) = §8.150E-01 
YERAC 9) = 1.000E+00 ;YFRAC 6) = 1.0Q0E+00 
YERAC 7) = 1.Q00E+00 ;YFRAC 8) = 8.200E-01 
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YFRAC 9) = 1.000E+00 ;YFRAC 10) = 1.190E+00 
YFRAC 11) = 1.000E+00 ;YFRAC IiZ). = 3.0508-01 
YFRAC 13) = 1.000E+00 ;YFRAC 14) = 1.520E+00 
YRRAC 15) = 1.0Q000E+00 ;YFRAC 16) = 3.050E-01 
YRRAC 17) = 1.000E+00 ;YFRAC 18) = 33050E-01 
YeRAC 19) = 1.000E+00 ;YFRAC 20) = 8.350E-01 
YFRAC 21) = 1.000E+00 ;YFRAC 22) = 4.600E-01 
YFRAC 23) = 1.000E+00 ;YFRAC 24) = 3.050E-01 
YFRAC ( 25) = 2.000E+00 ;YFRAC ( 26) = 2.780E+00 
YFRAC ( 27) = 1.0Q00E+00 ;YFRAC 28) = 3.660E+00 
YFRAC 274 =. .1.0008+00 ; c 30) = 1.000E+00 
KKK KKAKKAAKKKKKKRKKAKKKKKKAREKK KKK KHRKKKKKKKKKERKKAKAKKKKKK 


GROUP 5. Z-DIRECTION GRID SPACING 


PARAB = 9 
NZ = 3 
ZWLAST = 1.Q00E+00 

HETHOD OF PAIRS USED FOR GRID SETTING 
Be RaC lL) = -2 -OUCELUC Zr RAC 2) = IeeZ20E+00 
ZERAC Sy= 2.000E200 :ZFRAC = 2.440E+00 
GFRAC 5) = 2.000Et00 : ZFRAC = 1.780E+00 
ZFRAC 7) 5= |. OOOO > ZERAC = 3. 050E=01 
ZEFRAC 9) = 1.000E+00 ;ZFRAC 10) = 1.370E+00 
ZFERAC bee = 1 QOOUR BOO 2rRAC 12). =_33050E-01 
ZERAC 13) Parr OUCETO0 -ZrRAC 14) = 3.615E+00 
ZFRAC 15) = 1.000E+00 ;ZFRAC toy = 120002700 
ZERAC = “eee ee 00 > ZFRAC Te) = 3, 0S50E=0 
ZFRAC Poee— 1, 00UE2O0. -ZFRAC 20) = Pe220E.00 
ZFRAC 21) = 1.000E+00 ;ZFRAC 2235 =— 3,050E-01 
ZERAC 23)0= 1.0008700 -ZFRAC 24) = 3.680E+00 
ZE RAC 25) = 1.000E+00 ;ZFRAC 267) = 43570E5-01 
ZERAC 27) = 1.Q000E+00 ;ZFRAC Ze) = “S2000E+06 
ZERAC oe = 1.000E+00 ;ZFRAC 25 = 4.570E-01 
ZERAC = 1.000E+00 :ZFRAC = 3. 


een OOOE+00 
Be ECA TELIA ak rer Te rica EO Ria ia ade piceninoicininieiciek 
GROUP 6. BFC AND GEOHETRY CHANGES 


GROUP 7. VARIABLES: STORED,SOLVED,NAMED 
ONEPHS = a 

NAME( 1) =Pl ; NAME 
NAME( 5) =V1 ;NAME 
NAME(12) =KE ;NAME 
NAME(14) =H1 ;NAHE 
NAME(20) =NPOR ;NAME 
NAME(22) =VPOR ;NANE 
NAME(24) =TMP1l ;NAME 
SeLUPN(P1l , 
SOLUIM(Ul , 
SeLurNVvi , 
SUEomahl  , 
Sleek , 


=Ul 
=W1 
=EP 
=HPOR 
=EPOR 
=Viol 
=RHO1 


hb 
LO tl) ~J) ty 


NON PO 
Cito 


SOLUTN(EP 

SOLUTN(H1 

SOLUTN(HPOR 
SOLUTN(NPOR, 
SOLUTN(EPOR, 
SOLUTN(VPOR, 
SOLUTN(VIST, 
SOLUTN(TMP1, 
SOLUTN(RHO1, 
DEN1 = 
VIST 


EPOR ;HPOR = 19 ;NPOR = ZOmE VE OR = 22 
WEIL? 1 = 


4 
KK KAKI IK AAR AIK ARIK IRR R IR RIK IR IRI IR IR IR IIIA IR IRR 


Geeuree. "PERS &pDEVICES 
Plats (el ,Y,Y,Yie,¢%,N 


Veee eee ee eKKKe 
wee ee eee eee 
poreeeeeee eee 


BRO~ - ys es ye ye yy FMF 


WWI AQZABZAAGGs 


POM PO 
kK 


BEGHS (U1 /Y,¥,¥,N,Y,N 
Meeiowevl§=64y,Y,1%,N,Y,N 
tiara. ,Y,¢,1,N,7%,N 
TERHomCKkE ,Nyv;Y,N,Y,N 


By 


EP 
H1 


TERMS 
TERS 
Das 
GALA 
NEWRH1 
NEWENT 


} 


1.Q000E+00 


RAR RIA KAR RII KIN IR TE PPR TSE Te eo ee 
GROUP 9 PROPERTIES 


RHOL =-1.016E+04 ;TMP1l =-1 .013E704) emu =-1 .O15E+04 
RHOLA = 0.Q00E+00 ;RHOIB = 3.4962-03 7 RoOlC eee t-oL 
PRESSO = T.@ieErO5 

TMP1A = 1.000E-20 ;TMP1B = 9.960E-04 ;TMP1C = 0.000E+00 
TEMPO = Q.Q000E+00 

ELIA = 0, G00E.00 >, ELIB =.0'.000E+O0) ELIE = 0.Q000E+00 
ENUL = 1.000E-05 ;ENUT =-1.014E+04 

ENUTA = 0.Q000E+00 ;ENUTB = 0.Q000E+00 ;ENUTC = 0.000E+00 
DRHIDP =-LeGterree 

PRNDTL(U1 = 1.QQ00E+00 ;PRNDTL(V1 = 1.Q00E+00 
PRNDTL(W1 = 1.0Q00E+00 ;PRNDTL(KE = 1.000E+00 

PRNDTL( EP = 1.Q000E+00 ;PRNDTL(H1 = 1.000E+00 

PRE UL = 1.000ESGO SER V1 = 1.Q00E+00 

PAG W1 = 1.0Q0E+00 ;PRT KE = 1.Q00E+00 

Ena EE = 1.314E+0 Pa H1 = 1.000E+00 
KAKKKKKKARAKK AARAARL ARE RERKKKKKRKAKEKEKAKEKERKR RIK 


GROUP 10.INTER-PHASE TRANSFER PROCESSES 
KERKAKKAKKRAKKARKKKAAKRK AKA KKK RAKKAKKAAKRRKAKKKKERKAKARKERRRKKKRERAK 


GROUP 11.INITIALIZE VAR/POROSITY FIELDS 

PETC Pe I = -1.022E+04 ; NIT(U1 = -1.022E+04 
Pier Vv 1 = -],022ZE+O4 Bernie = -1.022E+04 
Pe CIE. = -1.022E+04 ;FIINIT(EP = -1.022E+04 
Be = -1.022E+04 | INIT (HeOR) =e OCU n we 
FIINIT(NPOR) = 1.QQ00E+tOG [RRENIT(EPOR) = el COGE reo 
FLIINIT(VPOR) = 1.Q00E+00 ;FIINIT(VIST) = ~eG2zetea 
FIINIT(TMP1) = =-1.022E+04 ;PLINTI(RHO]) 9= = fee7z2 ee 
KEYVELE = EB 

KELIN = 


0 
KRAEKARKRKRKKKRKEKRKRKRKKRKRERKARRKRKKRKRKRKRRKRKRRRKKRRRKRKRKKRRKRKERERERKKK 


GROUP 14. DOWNSTREAM PRESSURE FOR PARAB 
KRAKKKKKKKKKKKKAKKKRKKKKKKAKKKRKKRKEKKKKEKAKKKERKKKEKRKKKRKKKKKKKK 


GROUP 15. TERMINATE SNEEES 

LSWEEP = 550 = ise i 

LITHYD = 1 ;LITELA = ee he = 1 ;ITHCL = 
RESREP (EE = 1.Q00E-06 ;RESREF(U1 = 1.Q00E-06 
RESREF (V1 = 1.Q00E-06 ;RESREF(W1 = 1.000E-06 
RESREF (KE = J2O00-OGe Bron er (Er = 1.000E-06 
RESREP (ia = 


KRAAKKKKRAKK dekh Rick RRR AL AKI IRIR ARIA AAI EER 





GROUP 16. TERMINATE ITERATIONS 

Dien Cel = 20 ;LITER (Ul = ae 

Liner? (Vl =. 20" - fiver SS 

LITER (KE = 20 :LITER (EP = ee 

LITER (H1 =... 26 

ENDIT (Pl = 1.000E-06 ;:ENDIT (Ul = 1.000E-06 
ENDIT (Vl = 1.000E-06 :ENDIT (Wl = 1.000E-06 

a KE = 1.000E-06 ;ENDIT (EP = 1.000E-06 
ENDIT Hl od FOE TOO ee 
GROUP 17. RELAXATION 

RELAX(P1 ,FALSDT, 1-.000E-03 

RELAX(UlL ,FALSDT, 7.000E-02 

RELAX(V1 ,FALSDT, 7.000E-02 

RELAX(W1 ,FALSDT, 7.000E-02 

RELAX(KE ,FALSDT, 1.000E+10 

RELAX(EP ,FALSDT, 1.000E+10 

RELAX(H1 ,FALSDT, 1.000E+10 

RELAX(VIST,FALSDT, 1.000E+00 

RELAX(TMP1,FALSDT, 1.000E+00. 

RELAX(RHO1L.FALSDT. 1 OE+00 


CTO I TOIT IG EL eee 
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GROUP 18. LIMITS 

VARMAX(P1 = 2.500E+05 ;VARMIN 
VARMAX (U1 = 2.Q000E+02 ;VARMIN 
VARMAX (V1 = 5.000E+02 ;VARMIN 
VARHAX (W1 = 1.200E+03 ;VARMIN 
VARNAX (KE = 1.Q000E+06 ;VARMIN 
VARHAX( EP = 1.000E+06 ;VARMIN 
VARMAX(H1 = 2.008E+06 ;VARMIN 
VARMAX(HPOR) = 1.Q00E+10 ;VARMIN 
VARMAA(NPOR) = 1.Q000E+10 ;VARMIN 
VARMAA(EPOR) = 1.Q000E+10 ;VARMIN 
VARMNAX(VPOR) = 1.Q00E+10 ;VARMIN 
VARNAX(VIST) = 1.Q000E+10 ;VARMIN 
VARMAX(TMP1) = arias ; VARMIN 
VARMAX(RHO1) = 2.000E+00 - VARMIN 
a iene 
GROUP 19. EARTH CALLS TO GROUND 
USEGRD = T ;USEGRX = T 
NAMGRD =NONE 


Pl =-2.500E+05 
Ul =-2 .QQOE+02 
Vil =-5 .QO00E+02 
W1 =-2.000E+02 
KE = 1.000E-08 
EP = 1.000E-08 
Fn = 2.962E+05 
HPOR) =-1.Q000E+10 
NPOR) =-1.000E+10 
EPOR) =-1.000E+10 
VPOR) =-1.Q00E+10 
VIST) = 1.000E-20 
TMP1) = 2.950E+02 
RHO1) = 1.000E-Ol 
RERKRKKARRARRKERKRKREKRRRRKRKRRREEK 
STATION 


GENK = tf 
RRARKERAAKARRKAKARAKKKKAEKKAKAKKAARARAARRAAKARRAARAARARRARKKAKRARKKKRK 
GROUP 20. PRELININARY PRINTOUT 


ECHO 


a ae if 
RRERKRKARARKRRKRKRRRRRRRRKRKRRARRRRRRRRRREIAK 


GROUP 
INIFLD 
OUaE UT 
OUTPUT 
OUTPUT 
OUTPUT 
OUEEUT 
OUTPUT 
OUTPUT 
OUTPUT 
OUTPUT 
GUTPuUT 
OUnE UT 
OURrUT 
OUREUT 
GPTPur 


REERAKKRKRKRERARR 


GROUP 22. MONITOR PRiVT-OUT 
2 5 TYuer 


I AMON 
NPRMON 
HIGHLO 


21. PRINT-OUT OF VARIABLES 
= F ;:SUBWGR = F 
Pl ,Y.N,N,Y.¥,Y 

ie sy f,N,Y,Y,Y% 

Pe oF ,N,Y,Y,Y 

W1  {Y,/NUN,Y,¥,¥ 

Peer oN, YY YL 

EP ‘Y'N/N YY Y 

one N,N,Y,Y,Y 

PmeOR N,N,N,N,N,N 
NPOR,N,N,N,N,N,N 
EPOR,N,N,N,N,N,N 
VPOR,N,N,N,N,N,N 

heel, Y,uN,wWW,N 

iM), Y,N,N,N,N,N 

eee ee 


10000 :NPRMNT 


5 


ZMON 


el 
10000 : TSTSWP 


5 
1 


KREAKRKRRKAKAKAKARKRKAKKKRKRRRRRERERREA 


REKAKRKARAEKRKRERRRRRRRRRRRRRRERKRRRRRERERE 


RAH E RA RRA A AREA RAAKA RRR RARARRARRRARARAAAAKKE 


GROUP 
NPRINT 
NAPRIN 
bier all 
NZPRIN 
AZPR 
WeELIF 
ITABL 
ABSIZ 
NTZPRF 
Penk 


HOSP ALGHeoe Yel USED 
GROUP 24. BUIPS FOR RESTARTS 
SAVE 


NSAVE 
LUSAVE 


23. 


PIELD PRINT-OUT & PLOT CONTROL 
10000 ;NUNCLS = 5 
1 “Larne = te Lark = 6 
i ; IYPRE = c > 1YPRL = 16 
1 ;f2eke = i. -IZPRL = i 
EF 3Y2ER = E 
1 ;IPLTL = 250 ;NPLT = 15 
o77ir ROR 3 
5.000E-01 ;ORSIZ = 2.000E-01 
1 :NCOLPF = “50 
2 :NCOLCO = 45 ;NROWCO = 20 


T ;AUTOPS = 


=DFO9 


E 


FOR, THES GROUP 


;NOWIPE 


Ak GRID-GEOMETRY INFORMATION *** 
m-GGOnem@ieries OF THE CELL CENTRES 
Beeo2e +00 


DesOUEaO 1 


1. 205Er00 


4.0385E+00 
Y-COORDINATES OF THE CELL CENTRES 


TwOVE-02 5 


S/on=01 


1.465 


E+00 
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E 


2,9375500 


22ST SET O00 


3.432E+00 


Be oo0ET OU 


4.127E+00 5 .040E+00 Seve 2h Oo 64207 B00 6827200 
1 47 SELOO 1,8575200 9.400E+00 1. 216E 01 1.540E+01 
lady Seeor 
GZ-COORDINATES OF THE CELUSCE RES 
o- eek-O1 1.630E+Ge a, ooOEt00 6. L0GET00 S22Or +00 
JOO 1.1L0OsEtem dela (ls Onl 1227 Pee Od 1.467E+01 
eoo7 BOL 1. /6S35a00 ies Je TOE a eee 222Gb +01 
Ze 52 Leo 2.494E+01 Zaooi +01 2.840E+01 
--- INTEGRATION OF EQUATIONS BEGINS --- 
ie. oS Lee ee 1 SWEEP = 501 
TOTAL RESIDUAL/( 1.000E-06 ROR! Is 1.20s8708 
TOTAL RESIDUAL/( 1.Q000E-06) BOR UL is 4. 20a 02 
TOTAL RESIDUAL/( 1.0Q0B-GG maORT Vi iS 5.0ga 09 
TOTAL RESIDUAL/( FE.GO0E-06) FOR Swi TS Jest eed 
TOTAL RESIDUAL/( 1.00@E-06) S50Rere IS 2 6@36E412 
TOTAL RESIDUAL/( 1.00@0K-06) BaOReee IS 3.476E+14 
TOTAL RESIDUAL/( 1.000E 206 ) Wher eE IS 0.00@ETG0 
LIME “Smee. = t SWEEP = 502 
TOTAL RESIDUAL/( 1.000E-06) FOR PI IS 1.410508 
TOTAL RESIDUAL/( 1.000E-06) FOR U1 IS 5.263b 500 
TOTAL RESIDUAL/({ 1.000E-06) FOR V1 tae 5.910 se9 
TOTAL RESIDUAL/({ 1J000ER-C6) FORsnT IS. 2.22 7b 
TOTAL RESIDUAL/( I.QQ0E-06) FORGKE IS .3.500hrzZ 
TOTAL RESIDUAL/( 1 000E-06)} S2ORUEF IS 3.4928 714 
TOTAL RESIDUAL/( 1.000E-06) FOR H1 IS 3.244E+13 
TIME STEP = it SWEEP = 503 
TOTAL RESIDUAL/( 1.000E-06) FOR Pl IS l.S53b5e6 
TOTAL RESIDUAL/( 1.000E=06) Sher Ul IS 3.924E%02 
TOTAL RESIDUAL/( 1.000E-06) FOR V1 IS 4.l135Res 
TOTAL RESIDUAL/( 1.000E-06) FOR WI IS 1.4658 0n0 
TOTAL RESIDUAL/( 1.000E-06) FOR KE TS 3.0672 a 
TOTAL RESIDUAL/( 1 O00GE—-06) “FOR @EP IS 3.7638 ue 
TOTALVRESIDUAE){ 1000E=UG) Mier ens IS 3.24¢4Ra13 
LINES Stee = i SWEEP = 504 
TOTAL RESIDUAL/( 1. O0G0CE=06) Shona rl IS 1.6758202 
TOTAL RESEDUAL,( 1000-06 ) Baer aul TS. 5.352805 
TOTAL RESIDUAE,( 120005 -06 ) ShCRa 1 TS 5.93tenoe 
TOTAL RESIDUAL/( 1.000E-06) FOR Wl IS 2.223870 
TOTAL RESIDUAL/( 1.Q000E-06) FOR KE IS 1.87sbplz 
TOTAL RESIDUAL/( 1.000E-06) FOR EP IS. 3.38/72 
TOTAL RESIDUAL/( 1.000E-06) FOR H1 iS 2. Seems 
GE Slee = 1 SWEEP = U5 
TOTAL RESIDUAL/( 1.000E-06) FOR Pl IS less Enes 
TOTAL RESIDUAL/({ 1.000E-06) FOR U1 IS 4.270R 709 
TOTAL RESIDUAL/( 1.000E-06) FOR V1 IS 4.944E+09 
TOTAL RESIDUAL/( 1.000E-06) FOR W1 Foe 1 soSUbatO 
TOTAL RESIDUAL/( 1.000E-06) FOR KE Tom 1. 949E412 
TOTAL RESIDUAL/({ 1.Q000E-06) FOR EP Ise 3-S93Er14 
TOTAL RESIDUAL/( 1.000E-06) FOR H1 Poe 2.560 S5eglS 
EE Saer = Ih SWEEP = 506 
TOTAL RESIDUAL/( 1.Q00E-06) FOR Pl Is 1.0/;GRm0e 
TOTAL RESIDUAL/( 1.000E-06) FOR Ul IS -2:@528t03 
TOTAL RESIDUAL/( 1.Q00E-06) FOR V1 IS 4.394E+09 
TOTAL RESIDUAL/( 1.000E-06) FOR W1 15) 1 276E 10 
TOTAL RESIDUAL/( 1.000E-06) FOR KE [Sal oer te 
TOTAL RESIDUAL, 1.000R=06 | Shek ier IS 3 a@67Erl4 
TOTAL RESIDUAL/( 1.000E-06) FOR H1 IS 1.2a4B abs 
TAME STEP = 1 SWEEP) —s507 
TOTAL RESIDUAL/( 1.000E=06) FOR Pl IS 7. #ggn207 
TOTAL RESIDUAE/( 1.000206) FOR Ul IS 2.443E+09 
TOTAL RESIDUAL/( 1.000E-06) FOR V1 IS 22637obe09 
TOTAL RESIDUAL/( 1.000E-06) FOR Wl IS Pevlonri0 
TOTAL RESIDUAL/( 1.000E-06) FOR KE IS. LeZ2008+12 
TOTAL RESIDUAL/( 1.000E-06) FOR EP IS- 3.37@E+14 
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. 200E+00 
» 200E+00 
- 200E+00 
-200E+00 


MRRP RRR RRP RPP RHR RRP 


NNNNMNMMNNMMNNNNDNY 
in 


~.950E+02 
~.950E+02 
.950E+02 
-950E+02 
.950E+02 
.950E+02 
~950E+02 
~-950E+02 
-950E+02 
.950E+02 
~951E+02 
~951E+02 
oD LE+O2 
.951E+02 
.951E+02 
.951E+02 


.201E+00 
. 201E+00 
»201E+00 
~-201E+00 
»-201LE+00 
~-20LE+00 
roe +O0 
.201LE+00 
PoO LE +00 
.Z01LE+00 
AZOLE +00 
- 200E+00 
- 200E+00 
- 200E+00 
-200E+00 
- 200E+00 


TX= 


IY= 
IY= 
TY= 
IY= 
IY= 
i= 
1y= 
TYy= 
IY= 
IY= 
= 
i= 
IY= 
IY= 
IY= 
Iy= 
IX= 


— BRO WD 10) ~) 01 


-201LE+00 
.201E+00 
-201E+00 
-201E+00 
-2OLE+00 
-201E+00 
~.201LE+00 
-201E+00 
» 20LE+00 
,201E+00 
»>Z01LE+00 
-2Z01E+00 
~-201E+00 
-201E+00 
- 200E+00 
eee oo 


Pb Ph Pe Pe 


KRAEKKAAKAKKKKKKKKAKKAKKKKKKRKRKRKRKKKKRKKKKKRRKRKRKRKRKRKKRKRKKRKRRKRRERKRKRKKE 


TIME STP= 


FLOW FIELD AT ITHYD= 


FIELD VAEUBSe@e Fl 


= 16 -6.608E+00 
IY¥= 15 ~—( seer 00 
IY= 14 -8.487E+00 
IY= 13 > =Seeeeto0 
TY= 12 ‘=9eie 2700 
EY= 11° =S8ior+00 
TY= 10 =982266.,00 
Iy= 9 -=9.4eee+00 
TY¥= 8 =957s38700 
IY= 7 #955222 700 
TW= 6 ~s)2076E.0 1 
IY= 5 <-1.034E+01 
fy= 4 Fl OE TOl 
IY= 3 -8.840E+00 
iY¥= 2 eevee rO0 
i= 1 <2esoor-.00 
IXx= i 
TY= 16 =505/4E+00 
Ty¥= 15 -657052-00 
T¥= 14 #)].7626400 
T¥= 13 48.6405700 
T¥= 12 S38 7550.00 
Y= 11 =8.7225700 
TY= 10 -8.840E+00 
IyY= 9 -8.840E+00 
IY= 8 -8.836E+00 
i= 7 2.9m .00 
T= 6 2. soTeT00 
Ty¥= 5. =Ba70Z2e100 
TY¥= 4 -8.002E+0C 
Iv= 3 Gn eeet00 
IY= 2 =SegeeGR+00 
I¥= 1 =52se7E+00 
IX= 6 
BPIELD VALUBSSOr sl 
Iy= 16 “= eazlsE—02 
y= 15 =e -ol 
y= 14 © =e -01 
Iy= 13 =]02e=-01 
IY= 12 =Oaen-02 
Iv= ll =—Seison-02 
Iy= 10 =2 Ol 
X= 9 =Se30GE=-01 
IY= 8 -6.408E-01 
Iv 7  =6iecze-02 


1 SWEEP NO= 


550 ZSLAB NO= 2 ITERN NO= i 
1, IZ= 2, ISWEEP= 550, ISTEP= 
-6.60S5E+00 -6.423E+00 -6.257E+00 -6.082E+00 
-7.364E+00 -7.275E+00 -7.217E+t00 -7.161E+00 
-8.549E+00 -8.473E+00 -8.426E+00 -8.367E+00 
-9.481E+00 -9.421E+00 -9.389E+00 -9.332E+00 
-9.275E+00 -9.287E+00 <-9.341E+00 -9.368E+00 
-9.276E+00 -9.286E+00 -9.344E+00 -9.387E+00 
-9.322E+00 -9.311E+t00 -9.363E+00 -9.414E+00 
-9.419E+00 -9.349E+00 -9.347E+00 -9.419E+00 
-9.508E+00 -9.451E+00 -9.370E+00 -9.441E+00 
-9.741E+00 -9.903E+00 -9.802E+00 <-9.734E+00 
-1.074E+01 -1.080E+01 -9.888E+00 -9.436E+00 
-1.067E+01 -1.086E+0O1l -1.026E+O1l -9.771E+00 
-1.033E+01 -1L.1OOE+O1l <-1.029E+01 -9.371E+00 
-9.241E+00 -9.879E+00 -9.128E+00 -8.052E+00 
-6.816E+00 -8.101E+00 -7.325E+00 <-5.832E+00 
-6.096E+00 -7.955E+00 -7.246E+00 -5.534E+00 
2 5 4 5 
-2.700E-01 -3.967E-01 -5.002E-01 -5.098E-01 
-4.238E-01 -5.369E-01 -6.017E-01 -6.017E-01 
“4.384E-01 -5.699E-01 -6.569E-Ol1 -6.664E-01 
“4.252E-01 -5.778E-0Ol1 -6.914E-01 -7.107E-Ol 
“4.176E-01 -5.563E-01 <-6.580E-0O1 -6.817E-0O1 
-4.506E-01 -5.598E-01 -6.345E-0O1 -6.648E-01 
-5.778E-01 -5.853E-01 -5.891E-01 -6.336E-01 
-7.299E-01 -7.129E-01 -6.060E-01 -6.467E-01 
-8.280E-01 -7.983E-01 -6.703E-01 -6.757E-01 
-6.325E-0O1 -8.276E-01 -8.207E-0Ol1 -7.913E-01 


o2 


IY= 
IY= 
IY= 
Iv= 
IY= 
IY= 


LX= 


tH ho G),.AR WON 
Nr ANE FR 


ZU 28e0 | 
~-160E-01 
E277 e-OUl 
~3547E-01 
-355E+00 
aeo9e +00 


FELD VALUES OF V1 
T= 15 


Le 
IY= 
IYv= 
[Y= 
IYv= 
iyv= 
IYv= 
IY= 
LY= 
iv= 
LYy= 
ia 
LYv= 
IYv= 

PA 
IY= 
= 
i= 
TY= 
aE 
[Y= 
[v= 
[v= 
IY= 
[Y= 
IYv= 
TY= 
[v= 
ryr= 
IY= 


ITA= 


'etetttt 
NRRR~AIOF PRE NN WW ph 


MPMoW bOo~] COW 


-6.850E-01 
~515E-01 
.574E-01 
-650E-01 
~-129E-01 
-220E-01 
.943E-01 
~132E-01 
.767E-01 
Pole +00 
-626E+00 
-070E+00 
.O73E+00 
.554E+00 
~LO8E-01 


i 

.877E-0O1l 
.979E-01 
.349E-0Ol 
-461E-01 
wees 0 1 
444E-01 
~557E-02 
~LO2E-02 
.O89E-01 
-426E-01 
~639E-01 
~-237E+00 
-400E+00 
.L62E+00 
.O46E-01 


FIELD VALUES OF W1 
IY= 16 


TY= 
lY= 
lw ge 
Iv= 
[v= 
[Y= 
IYv= 
IYv= 
Iv= 
i= 
[Y= 
[Y= 
IYv= 
IY= 
Iv= 


i 
rNWAROoO~)OWO 


PRR RRR E 
WANnIDOOFRNWAMNOA 
t F ~) 01 01 01 01 OF 01 OV ni Ui) 

re 


5.766E-01 
8.881E-01l 
.074E+00 
. 166E+00 
-143E+00 
.102E+00 
.037E+00 
~119E+00 
-186E+00 
-442E+00 
.635E+00 
.793E+00 
. 704E+00 
-630E+00 
2-016E+00 
~-524E+00 


.693E-01 
~397E-01 
~478E-01 
~994E-01 
~-529E-01 
.342E-01 
.2/73E-0l 
6326-01 
.888E-01 
27 ee -01 
.057E+00 
~124E+00 


ce I ll ere en 


PSO 9e=01 
.045E-01 
.570E-01 
-200E+00 
.183E+00 
.894E+00 


-064E-01 
-201E-O1 
Ole 01 
»412E-01 
oor = Ol 
ween Oi) 
soguE Ol 
su gon- Oil 
-ogg2=01 
oe uO 1 
eleze =O) 
erOr=O I 
eee 0) 
oo B=O] 
same 


Jooue=—0] 
.284E-01 
.134E+00 
eZee +O0 
~-214E+00 
-FS0E+00 
~124E+00 
.189E+00 
'Zo0LT00 
-2Z10E-01 
~433E-01 


Sooo e-O2 


=o aoe OZ 
=a oon Ol 
-1.461E+00 
alle 


DHMH OAOAWH-H 
WW 


PO PRRPRPRPRPRHPRPHPEPPRPWOUM 
UJ 


rel elE = (© 1 
.230E-01 
.860E-01 
-030E+00 
-071E+00 
eee OO 


pao vee Ol 
ecoce ~Ol 
aoe =O] 
poe =O1 
.040E-01 
FoerE-01 
727 0E-04 
ome =O) 
sewer =O 1 
soeoe O01 
poouUE-O1 
-784E-01 
oom -O 1 
soo 92-01 
~o oe =02 


./90E-01 
~L97E-0O1 
-L36E+00 
~Z251E+00 
-262E+00 
~252E+00 
-240E+00 
~275E+00 
-294E+00 
-2387E+00 
-480E+00 
-456E+00 
-226E+00 
-550E-01 
-LOOE-01 
7425700 


33 


Mr 
if 


'oeoeesea 
po) OD oe ee 


CO PRR HRP RRRR REE OM 
if 


»265E-01 
-277E-O1l 
~992E-01 
.656E-01 
.438E+00 
.908E+00 


.967E-01 
-O10E-01 
-256E-01 
.570E-01 
~LOOE-O1 
.880E-01 
-524E-02 
.603E-01 
~546E-01 
-165E+00 
.269E+00 
-399E+00 
- 143E+00 
.664E-01 
298E-01 


.399E-01 
-419E-01 
-038E+00 
- L42E+00 
. L44E+00 
-124E+00 
~L22E+00 
.186E+00 
-2260E+00 
-285E+00 
-629E+00 
.677E+00 
.634E+00 
-003E+00 
~O81lE-0O1l 
~451E-01 


oe 


OME OU OOOOWOW OP 
(in 


NRPRRrOrR A 
in 


mo7Se=-01 
-405E-01l 
~553E-01 
-855E-02 
.014E+00 
page eT 00 


.906E-01 
~177E-01 
-443E-01 
-214E-01 
-776E-01 
~912E-01 
-306E-02 
-870E-02 
.573E-01 
.805E-01 
.078E+00 
~-575E+00 
-612E+00 
eZ OE+00 
-326E-01 


~564E-01 
peel 
me7 Le =-O1 
~905E-01 
.563E-01 
.342E-01 
.302E-01 
.733E-01 
-O05E-01 
-S29E-01 
-343E+00 
-407E+00 
~-487E+00 
.097E+00 
~754E-01 
-461E-01 


Iy= 4 1.222E+00 
iy=- .3 1.094E+00 
y= Zz 5.814E-01 
Iyv= 1 3.795E-01 
IXk= 6 
FIELD VALUES OF KE 
IY= 16 3.014E+00 
Iy= 15 1.081E+00 
IY= 14 1.027E+00 
IY= 13 1.074E+00 
Py= 12 9.180E-01 
Py=" 11 4.809E-01 
IY= 10 1.964E-01 
Iy= 9 1.515E-01 
Py¥= 6 1.411E-01 
Ivy= 7 1.409E+00 
IYy= 6 1.552E+00 
Iy= 5 1.967E+00 
Iy= 4 2.189E+00 
IY= 3 2.371E+00 
Iv= 2 3.593E+00 
Ivy= 1 9.283E-02 
IX= il 
IY= 16 4 .013E+00 
Ivy= 15 3.089E+00 
Ivy= 14 3.343E+00 
hye 3 3.571E+00 
Ivy= 12 3.468E+00 
IY= 11 3.437E+00 
IY= 10 3.421E+00 
Iy= 9 3.418E+00 
Ty= Gés 3.418E+00 
ly= @ 3.428E+00 
Iv= 6 3.291E+00 
Iv= 5 3.232E+00 
Iy= 4 2.928E+00 
Iy= 3 2.654E+00 
IY= 2 225725100 
ye il Seno 
FIELD VALUES OF EP 
IY= 16 5 .473E+00 
y= "15 1.250E+00 
IY= 14 1.180E+00 
Ey= 13 1.248E+00 
IY= 12 9.038E-01 
IY= ll 3.726E-01l 
Iv= 10 1.154E-01 
Iy= 9 7.985E-02 
Iy= 8 TeaZa5-02 
IY= 7 1.268E+00 
IY= 6 1.414E+00 
Iy= 5 1.737E+00 
Iy= 4 1.873E+00 
iy= 3 1.910E+00 
Iy= 2 S225 540 0 
Iy= 1 1.424E-01 
IX= iL 
Iv= 16 7.764E+00 
ITYy= 15 5 .604E+00 
Ivy= 14 6.179E+00 
IY= 13 6.674E+00 
i7= 12 6 .512E+00 
Iv= 11 6 .470E+00 
Ivy= 10 6.462E+00 
Iy= 9 6.459E+00 
IYy= 8 6.459E+00 
IYy= 7 6.502E+00 
IY= 6 6.269E+00 
Iy= 5 6.196E+00 


CONF OMNINRRNMORRHR W 
NO 


MPrROPRWHRRPRPRHWORRPRU 
NO 


20 Ps coo 
.120E+00 
-080E+00 
Soom EoO 
SP SUS Sr 
wee — OL 
woo Ol 
ooo L 
roa =O 1 
souve-O1 
soon =O 
poose +O 1 
oe eu 
-463E +00 
.384E+00 
eo e-UZ 


-469E+00 
.298E+00 
-244E+00 
-324E+00 
-590E-01 
.934E-01 
~-450E-01 
~-131E-01 
~112E-01 
~-462E-01 
~392E-01 
~O87E-01 
.694E-01 
polZe-Ol1 
solZzs+OO 
-225E-01 


2.978E+00 
1.215E+00 
Pez 5e+O0 
1.306E+00 
1.124E+00 
6.293E-01 
4.024E-01 
4.033E-01 
4.296E-01 


6.131E-01 
eo JE tOU 
ZeeZon too 
Th 02 


539 5E+00 
1.385E+00 
1.373E+00 
1.493E+00 
1.094E+00 
~O34E-01 
.982E-01 
~117E-01 
-469E-01 
.622E-01 
-695E-01 
-442E-01 
.861E-01 
.323E-01 
.423E+00 
-630E-02 


WOFRMOANNMNWWWNHWM 
W 


54 


3.028E+00 
1.425E+00 
1.486E+00 
io lPe+O0 
1.398E+00 
1.003E+00 
8.237E-01 


Ta oy BOL 
Jae ere 0 1 
1.324E+00 
1 ose reo 
oO 


5 #895E700 
1 J643E+00 
1.718E+00 
leo oZe 00 
1.482E+00 
9. /8e5-01 


1.214E+00 
as =O2 


NRRRNNWNNNNWWWNHUI 
oa 


3.318E+00 
2.054E+00 
2.225E+00 
Z2-407E+00 
Ze2ole+o0 
Zee ood 
2-J81E+00 
Ze 
2 
Z 
1 
i 
1 
A 
1 
5 


094E+00 


-LZ1E+00 
~191E+00 
.940E+00 
.815E+00 
.404E+00 
-480E+00 
-905E+00 
peewee 


~912E+00 
.818E+00 
~125E+00 
-441E+00 
.110E+00 
.974E+00 
.871E+00 
.902E+00 
.964E+00 
-153E+00 
.816E+00 
.607E+00 
-692E+00 
-LOZE+00 
-260E+00 
-899E-01 


Iv= 4 5.611E+00 
Ig= 3 4.937E+00 
IvY= 2 4.379E+00 
1 I ee 
FIELD VALUES OF H1 
a= 16 2.962E+05 
fv 15 2.962E+05 
IY= 14 2.962E+05 
IY= 13 2.962E+05 
T= I2 2.962E+05 
IY= 1l 2.962E+05 
IY= 10 2.962E+05 
iv 9 2.962E+05 
IYy= 8 2.962E+05 
i= 7 2.963E+05 
TY= 6 2.963E+05 
i= 5 2.964E+05 
iv 4 2.964E+05 
IyY= 3 2.963E+05 
i= 2 2.963E+05 
I¥= 1 2.963E+05 
IX= 1 
T¥= 16 2.962E+05 
y= 15 2.962E+05 
IY= 14 2.962E+05 
IY= 13 2.962E+05 
re 12 2.962E+05 
IY= ll 2.962E+05 
Iy= 10 2.962E+05 
TY= 9 2.96ZE+05 
im 68 2.962E+05 
ir—- 7 2.962E+05 
IY= 6 2.962E+05 
ee 5 2.962E+05 
Ivy= 4 2.963E+05 
IY= 3 2.963E+05 
a= 2 2.963E+05 
Ty= 1 2.963E+05 
IX= 6 
FIELD VALUES OF TMP1l 
IY= 16 2.950E+02 
TvY= 15 2.950E+02 
IY= 14 2.950E+02 
IY= 13 2.950E+02 
IY= 12 2.950E+02 
i= 1 2.950E+02 
Iv= 10 2.950E+02 
I¥= 9 2.950E+02 
iS 6 2.950E+02 
IY= 7 2.951E+02 
IY= 6 2.952E+02 
i= 5 2.952E+02 
IY= 4 2.952E+02 
Iy= 3 2.951E+02 
Iv= 2 2.951E+02 
iz 1 2.951E+02 
IX= 1 
IY= 16 2.950E+02 
Iig= 15 2.950E+02 
Ivy= 14 2.950E+02 
rye 13 2.950E+02 
iw 12 2.950E+02 
IY= ll 2.950E+02 
IY= 10 2.950E+02 
Iv= 9 2.950E+02 
Iv= 8 2.950E+02 
Lra= 7 2.950E+02 
Iv 6 2.950E+02 
Iv= 5 2.950E+02 


NNNMNMNMNNNNNNMNNMNNMND 
N) 


NMNMNNMNMMNMNMNM NPN A 
NO 


-962E+05 
~-962E+05 
-962E+05 
~-962E+05 
-962E+05 
~962E+05 
.962E+05 
~.962E+05 
~962E+05 
~963E+05 
.963E+05 
-963E+05 
.963E+05 
.963E+05 
.963E+05 
.963E+05 


.950E+02 
.950E+02 
-950E+02 
~950E+02 
-950E+02 
~950E+02 
~-950E+02 
~950E+02 
-950E+02 
Pool eroZ 
~951E+02 
eee FOZ 
951 E+02 
~951E+02 
~9I51E+02 
951 E+02 


MPM NMNMNNNNMPNNNMNMMNA) 
WJ 


NNMNMNMNMNMNMNMDNNDNNMNMNMNMN NN 
ty 


~-962E+05 
-962E+05 
-962E+05 
-962E+05 
-962E+05 
.962E+05 
~-962E+05 
~962E+05 
~963E+05 
-963E+05 
~963E+05 
-963E+05 
-963E+05 
.963E+05 
.963E+05 
~-963E+05 


~.950E+02 
-950E+02 
.950E+02 
-950E+02 
-950E+02 
.950E+02 
-950E+02 
~951E+02 
~.951E+02 
~951E+02 
-951E+02 
.951E+02 
~951E+02 
951E+02 
.951E+02 
~951E+02 


a0 


pr eaN NS De oa eee 


FOR Oneness Ce 


.962E+05 
P96ZE+05 
.962E+05 
-962E+05 
.962E+05 
.962E+05 
-962E+05 
~-962E+05 
-962E+05 
~962E+05 
.963E+05 
.963E+05 
-963E+05 
-963E+05 
.963E+05 
.963E+05 


.950E+02 
.950E+02 
.950E+02 
.950E+02 
.950E+02 
~950E+02 
.950E+02 
~950E+02 
951 E+02 
~951LE+02 
~951E+02 
-951E+02 
~95LE+02 
-951E+02 
~951E+02 
~951E+02 


NMMNMNNNMNNMNNNMDN MND Nh 
um 


NDMMNNM NN DNDMDN NPN 
Nn 


.962E+05 
.962E+05 
~962E+05 
.962E+05 
.962E+05 
.962E+05 
.962E+05 
-962E+05 
.962E+05 
.962E+05 
.963E+05 
.963E+05 
.963E+05 
.963E+05 
-963E+05 
-963E+05 


.950E+02 
.950E+02 
.950E+02 
-950E+02 
.950E+02 
.950E+02 
.950E+02 
-950E+02 
-.950E+02 
.950E+02 
~.9SLE+02 
~951E+02 
~951E+02 
~951E+02 
~951E+02 
-951E+02 


Tv= 4 2.951E+02 
TY= 3 2.951E+02 
ty= 2 2.951E+02 
Iv= 1 2.951E+02 
LxX= 6 
FIELD VALUES OF RHO1 
IY= 16 1.Z201E+00 1.201E+00 1.201E+00 1.Z201E+00 1.201E+00 
hy= 15 1.201E+00 1.201E+00 eZ OLe+O0 1.201E+00 1.201E+00 
IY= 14 1.201E+00 1.201E+00 1.201E+00 1.Z201E+00 1.201E+00 
y= 13 1.201E+00 1.201E+006 1.201E+00 1.201E+00 1.201E+00 
Lry= WZ 1.201E+00 1.Z01LE+00 1.201E+00 1.201E+00 1.201E+00 
IY= ll 1.201LE+00 eZore +00 1.2Z01LE+00 face +O0 1.2Z01E+00 
IY= 10 1.201E+00 1.201E+00 1.201E+00 1.201E+00 1.Z201E+00 
IYy= 9 1.201E+00 1.201E+00 IzeZOLE+O0 12 201E+00 1.201E+00 
Iy= 8 1.201E+00 1.201E+00 1.201E+00 1.201E+00 1.201E+00 
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we 
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Jee 
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201E+00 


-2Z01E+00 
~-2Z01LE+00 
-201E+00 
-200E+00 
- 200E+00 
- 200E+00 
- 200E+00 
oe 


Iz=— 13 ioe +00 
TY= 12 1.201E+00 
IY= 11 1.201E+00 
1¥= 10 1.201E+00 
Iy= 9 1.201E+00 
im= © 1.201E+00 
igs= 7 1.201E+00 
IY= 6 1.Z01E+00 
i= 5 1.201E+00 
IY= 4 1.201E+00 
Ty¥= 3. 1.200E+00 
IY= 2 1.200E+00 
m= WeZ00E+00 
IX= 
KRAERAKRAKKRKKKKKEKKKKRKKKK KAKA KKKKKKRK RK KKK KKK KKK KKKKKKKKKK 
TINE STP= 1 SWEEP NO= 550 ZSLAB NO= 5 ITERN NO= Ih 
FLOW FIELD AT ITHYD= 1, IZ= 5, ISWEEP= 550, ISTEP= 
FIELD VALUES OF Pl 
IY= 16 -8.620E+00 -8.642E+00 -8.460E+00 -8.276E+00 -8.049E+00 
TY= 15 -9.359E+00 -9.487E+00 -9.536E+00 -9.582E+00 +-9.591E+00 
ITY= 14 -1.444E+01 -1.461E+O1 -1.444E+0i -1.435E+01 -1.422E+01 
IY= 13 -2.404E+01 -2.426E+01 -2.351E+0O1 -2.303E+O1 -2.249E+0O1 
ITY= 12 -2.148E+01 -<2.297E+01 -2.500E+01 -2.684E+01 -2.739E+01 
IY= 11 -8.911E+00 -1.096E+01 -1.540E+01 -2.070E+0O1 -2.951E+0O1 
IY= 10 9.249E+00 4.812E+00 -5.428E+00 -1.718E+0O1 -3.196E+0O1 
Ivy= 9 6.724E+00 7.288E-01 <-1.079E+O01 -2.228E+0O1 -3.494E+01 
IY= 8 3.405E+00 -3.197E+00 -1.481E+01 -2.631E+0O1 -3.458E+01 
i= 7 1.129E+01 1.046E+00 -1.311E+01 -2.447E+01 -4.111E+01 
ITY= 6 -5.512E+01 -6.392E+01 -6.232E+01 -5.712E+0O1 -4.735E+01 
ITY= 5 -7.141E+01 -8.163E+0O1l -7.852E+01 -6.750E+01 -5.929E+01 
TY= 4 -8.015E+01 -8.539E+O1 -7.763E+O1l -6.756E+01 -6.177E+01 
IY= 3 -7.809E+O01 -8.343E+01 -8.020E+01 -6.877E+01 -6.143E+01 
IY= 2 -5.101E+01 -6.357E+0O1 -7.494E+01 -6.326E+01 -5.483E+01 
IY= 1 -1.439E+01 -4.965E+01 -7.254E+01 -6.212E+01 -5.368E+01 
IX= a Zz S 4 5 
IY= 16 -7.284E+00 
IY= 15 -8.942E+00 
ITY= 14 -1.279E+01 
Myv= 12 <=-1).913E+01 
Mm="12 -<2.307E+01 
IYV= 11 2-2.447E+01 
fe= 1@ =-2.5677E+01 
IY= 9 -2.863E+01 
IY= 8 -2.957E+01 
IY¥= 7 -3.317E+01 
IY= 6 -3.664E+01 
ITY= 5 -4.027E+01 
IY= 4 -4.163E+01 
TY= 3 -4.133E+01 
IY= 2 -3.917E+01 
IY= | “ae Od 
FIELD VALUES OF Ul 
IY= 16 -4.983E-02 -2.975E-O1 -4.681E-01 -6.092E-01 -6.264E-01 
IY= 15 -1.293E-01 -4.802E-01 -6.027E-01 -6.880E-01 <-7.085E-0O1 
IY= 14 -1.895E-01 -6.669E-01 -8.345E-01 -9.613E-01 -9.968E-0O1 
IY= 13 -2.227E-01 -8.710E-O1 -1.144E+00 -1.382E+00 -1.453E+00 
Iv= 12 -7.274E-02 -6.320E-O1 -7.627E-0O1 -1.146E+00 -1.476E+00 
IY= 11 -4.721E-01 -1.819E+0O0 -1.881E+00 -1.132E+00 -1.518E+00 
IY= 10 -9.765E-01 -2.674E+00 -2.280E+00 -5.634E-01 -1.309E+00 
IY= 9 -1.072E+00 -2.886E+00 -2.462E+00 -7.403E-01 -1.355E+00 
IvY= 8 -1.061E+00 -2.688E+00 -1.655E+00 1.748E-03 <-1.130E+00 
IY= 7 -1.443E+00 -3.570E+00 -2.859E+00 -5.867E-01 -1.454E+00 
IY= 6 -1.501E+00 -3.737E+00 -2.885E+00 <-1.669E+00 -1.884E+00 
IvY= 5 -1.249E+00 -3.735E+00 -3.725E+00 -3.528E+00 -3.420E+00 
IvY= 4 -1.251E+00 -3.677E+00 -3.971E+00 -3.867E+00 -3.763E+00 
IY= 3 -6.828E-01 -2.419E+00 -3.262E+00 -3.579E+00 -3.533E+00 
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IY= 2 1.235E+00 
LYy=) 1 3.901E+00 
IX= i 
FIELD VALUBS OF V1 
TY¥= 15 =885390E-01 
IY= ie -2.349E+00 
TY= 13 -3.191E+00 
IY= 12 -1.470E+00 
IY= 11 <-7.608E-01 
IY= 10 -2.837E+00 
IY= 9 2-5.493E+00 
IY= 8 -+6.429E+00 
IY= 7 -6.859E+00 
IY= 6 -1.163E+01 
IY= 5 -1.246E+01 
lTY= 4 2-1.247E+01 
ITY= 3 -1.145E+01 
IY= 2 -=8.200E+00 
LTY= 1 2*2.001E+00 
IX= 1 
T¥= 15 =620675-01 
ITY= 14 -9.439E-01 
TY= 13 -1.124E+00 
TY= 12 -9.570E-01 
ITY= 11 -9.905E-01 
ITY= 10 -9.399E-01 
IY= 9 -8.729E-01 
ITY= 8 -8.402E-01 
ITY= 7 -8.777E-01 
IY= 6 7-8.829E-01 
TY= 5 -+9.322E-01 
ITY= 4 -1.899E-01 
IY= 3 3.469E-01 
IY= 2 7.283E-01 
y= if 1.484E-01 
IX= 6 
FIELD VALUES OF W1 
Y= 6G 3.457E-01 
IY= 18 dee BoE ROG 
IY= 14 2.476E+00 
TY= 13 5.901E+00 
IY= 12 8.390E+00 
IY= 11 -6.499E+00 
ITY= 10 -6.499E+00 
IY= 9 -6.500E+00 
ITY= 8 -3.596E+00 
ITY= 7 -6.499E+00 
IY= 6 -+4.916E+00 
IY= 5 1.844E+00 
IY= 4 6.12Z21E+00 
IYy= 3 1.017E+01 
Iv= 2 1.294E+01 
Iv= 1 1.420E+01 
IX= il 
IY= 16 1.531E-01 
ty= 15 3.187E-01 
IY= 14 8.783E-01 
Ty= 13 2.043E+00 
hy=*) 2 2.736E+00 
IY= 1l 3.544E+00 
IY= 10 4.608E+00 
Iy= 9 5 .034E+00 
Ivy= 8 5 .251E+00 
TY= 7 5.803E+00 
IY= 6 5.412E+00 
Ty¥= & 4.872E+00 
Iy= 4 4.621E+00 
IY= 3 4.600E+00 
IY= 2 4.809E+00 
Ivy= 1 4.813E+00 
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FIELD VALUES OF KE 
IY= 16 
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IY= 4 
m= 3 
i= 2 
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~350E+01 
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.873E+00 
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-321E+00 
~955E+00 
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1.367E+01 
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FIELD VALUES OF Hl 


TY=76 2.962E+05 
T¥= 25 2.962E+05 
IYy= 14 2.962E+05 
TY= 13 2.962E+05 
ITY= 12 2.963E+05 
LY= ai 2.968E+05 
LY= 20 Z2.969E+05 
IYy= 9 2.969E+05 
Iy= 8 2.969E+05 
TYy= 7 2.968E+05 
IY= 6 2.968E+05 
Ty= 5 2.967E+05 
Iy= 4 2.966E+05 
ry= 92 2.964E+05 
Iv= 2 2.963E+05 
y= <1} 2.963E+05 
IX= 1 
IY= 16 2.962E+05 
tyv= 15 2.962E+05 
IY= 14 2.962E+05 
TY= 13 2.962E+05 
ITY= 12 2.962E+05 
ITY= 11 2.962E+05 
TY= 10 2.962E+05 
IYy= 9 2.962E+05 
IYy= 8 2.962E+05 
IY= 7 2-9625+05 
Iy= 6 2.962E+05 
LY= 5 2.962E+05 
Iy= 4 2.962E+05 
Ly= =3 2.962E+05 
Iv= 2 2.963E+05 
ee 1 2.963E+05 
x= 
FIELD VALUES OF TMP1l 
ITY= 16 Z2e950E +02 
rY= 15 2.950E+02 
IY= 14 2.950E+02 
TY= 13 2.950E+02 
IY= 12 2.951E+02 
IY= 1l 2.956E+02 
ITY= 10 2.957E+02 
Iy= 9 2.957E+02 
IY= 8 2.957E+02 
IY= 7 2.956E+02 
IvY= 6 2.956E+02 
b= 2.955E+02 
Ivy= 4 2.954E+02 
Ivy= 3 2.952E+02 
IY= 2 2.951E+02 
hy= “] 2.951E+02 
LX= 1 
IY= 16 2.950E+02 
IvY= 15 2.950E+02 
IY= 14 2.950E+02 
Ivy= 13 2.950E+02 
IY= 12 2.950E+02 
IY= 1l 2.950E+02 
IY= 10 2.950E+02 
IYy= 9 2.950E+02 
Iv= 8 2.950E+02 
IY= 7 2.950E+02 
IY= 6 2.950E+02 
Iy= 5 2.950E+02 
ITy= 4 2.950E+02 
IY= 3 2.950E+02 
IY= 2 2.951E+02 
Iy= 1 2.951E+02 
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FEEED VALUES OF RHO1 


IY= 16 1.20 2E+00 1.201E+00 1.Z01E+00 1.201E+00 1.201E+00 
Iy¥= 15 1.201E+00 1.201E+00 1.2Z01E+00 eZ eae) e220 le +O0 
IY= 14 eZee +O0 1.201E+00 1.201E+00 1.201E+00 eZ OLE+O0 
te 13 1.Z201E+00 1.201E+00 1.201E+00 1.201E+00 ieee be +00 
IY= 12 eZ00E+00 1.200E+00 1.201E+00 1.201E+00 20 le OO 
lY= ll 1.198E+00 1.199E+00 1.200E+00 ior E+OO ieee bE Oo 
IY= 10 1.198E+00 1.Z200E+00 e200 E +00 1.Z201E+00 eZ 00E+O0 
f= 9g 1.198E+00 200700 1.200E+00 1.200E+00 1.Z00E+00 
IYy= 8 1.198E+00 1.200E+00 1.200E+00 1.200E+00 1.200E+00 
Iy= 7 1.198E+00 i Z00E+00 1.200E+00 1.200E+00 1.200E+00 
IY= 6 1.198E+00 1.199E+00 1.200E+00 ZOO +O0 1.200E+00 
i= 5 1.198E+00 1.199E+00 1.200E+00 PaZ0OE+00 1 .Z200E+00 
Iy= 4 1.198E+00 1.199E+00 1.200E+00 1.200E+00 1.200E+00 
Iy= 3 1.199E+00 1.200E+00 1.200E+00 1.200E+00 1.200E+00 
k= 2 1.200E+00 1 .200E+00 1.200E+00 1.200E+00 1.200E+00 
il “ae 1.200E+00 meee OO lee CoE +o 1.200E+00 
Iy= 16 1.201E+00 
i 15 1.201E+00 
IY= 14 Pez +00 
Ey=— 13 1.201E+00 
iz= £2 1.20 LE700 
IY= ll 1.201E+00 
IY= 10 1.201E+00 
m= 9 1.201E+00 
Iy= 8 1.Z01LE+00 
IyY= 7 1.200E+00 
IY= 6 1.200E+00 
y= 5 1.200E+00 
Iv= 4 1.Z00E+00 
Iv 3 1.200E+00 
IY= 2 1.200E+00 
IYy= 1 1.200E+00 
IX= 6 
KAKA KK IK IRKAARIKI IK IKIKRARKAKAKAAAKKRKREKEKARAKRAKR AK A AKA AKA K 
TIME STP= 1 SWEEP NO= 550 ZSLAB NO= 6 ITERN NoO= 1 
FLOW FIELD AT ITHYD= 1, IZ= 6, ISWEEP= 550, ISTEP= 1 
PIED VALUES OF Pl 
IY= 16 -9.132E+00 -9.147E+00 -8.957E+00 -8.768E+00 -8.526E+00 
IY= 15 -9.467E+00 -9.676E+00 -9.891E+00 -1.002E+01 -1.011E+01 
IY= 14 -1.790E+O01 -1.810E+O1 -1.779E+0O1 -1.759E+01 -1.735E+01 
IY= 13 -4.656E+O1l -4.691E+01 -4.243E+01 -3.949E+01 -3.709E+01 
IY= 12 -+1.627E+02 -1.606E+t02 -1.284E+02 -1.040E+02 -7.788E+01 
IY= 11 -1.875E+02 -1.874E+02 -1.526E+02 -1.310E+02 -9.693E+01 
IY¥= 10 -1.643E+02 <-1.730E+02 -1.542E+02 -1.515E+02 -1.257E+02 
IY= 9 -1.498E+02 -1.648E+02 -1.557E+02 -1.602E+02 -1.393E+02 
IY= 8 -1.461E+02 -1.647E+02 -1.606E+t02 -1.655E+02 -1.475E+02 
Iv= 7 -+1.499E+02 -+1.826E+02 -1.910E+02 -2.066E+02 -1.802E+02 
IY= 6 2*1.531E+0O2 -1.837E+t02 -1.923E+02 -1.863E+02 -1.634E+02 
Iv= 5 -1.656E+02 -1.939E+t02 -1.845E+02 -1.653E+02 -1.502E+02 
IY= 4 -1.876E+02 -2.034Et02 -1.791E+0O2 -1.575E+02 -1.444E+02 
ITY= 3 -2.095E+02 -2.072Et02 -1.764E+02 -1.535E+02 -1.398E+CZ 
TY= 2 -1.952E+02 -1.959E+02 -1.655E+02 -1.458E+02 -1.316E+02 
Iv= 1 -1.724E+02 -1.901Et02 -1.637E+02 -1.451E+02 -1.306E+02 
Ix= iL Z e 4 5 
IY= 16 -7.681E+00 
TY= 15 2+9.472E+00 
TY= 14 -1.528E+01 
W= 13 =2.-625E+01 
TY= 12 -4.421E+01 
Iy= 11 =5.2047E+01 
IvY= 10 -6.100E+01 
IvY= 9 -6.689E+01 
IvY= 8 -7.026E+01 
IY= 7 -8.093E+01 
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IY= 6 -8.Q00E+01 
IY= 5 -7.936E+01 
T¥= 4 7 3855E cel 
IY= 3 2-7.729E+01 
IY= 2 -7.507E+01 
IY= 1 -7.502E+01 
IX= 6 
FIELD VALUES OF U1 
IY= 16 -5.277E-02 
IY= 15 -6.874E-02 
IY= 14 <-1.635E-01 
IY= 13 -4.486E-01 
IY= 12 -5.247E-01 
IY= 11 -8.230E-01 
IY= 10 -1.641E+00 
I¥= 9 -=1.930E+00 
IY= 8 -2.231E+00 
IY= 7 2-2.793E+00 
IY= 6 -3.072E+00 
IY= 5 -2.859E+00 
IY= 4 -3.066E+00 
TY¥= 3 -=3.205Et00 
IY= 2 -2.369E-01 
Iy= 1 a 
FIELD VALUES OF V1 
IY= 15 -8.906E-0O1 
IY= 14 -3.254E+00 
IY= 13 -5.429E+00 
IY= 12 <-1.106E+01 
Ev= ll =le2beeeo.t 
IY= 10 <-1.061E+01 
IY= 9 -9.551E+00 
IY= 8 -9.335E+00 
IY= 7 -+9.776E+00 
TY= 6 -8.774E+00 
ITY= 5 -8.730E+00 
IY= 4 -8.832E+00 
TY= 3 -8.885E+00 
ITY= 2 -6.466E+00 
IY= 1 -1.661E+00 
IXk= 1 
IY= 15 -6.307E-0Ol 
IY= 14 -1.247E+00 
IY= 13 -1.726E+00 
IY= 12 2-2.424E+00 
ITY= 11 -=2.535E+00 
T¥= 10 =2.320E+4+00 
IY= 9 -1.662E+00 
IY= 8 -1.468E+00 
IY= 7 -1.312E+00 
IY= 6 -2.691E-01 
IY= 5 -1.892E-01 
Iy= 4 6.413E-02 
Ivy= 3 2.814E-01 
Iv= 2 4.470E-01 
Iv= 1 aes 
FIELD VALUES OF Wl 
Iv= 16 3.024E-01 
he— 5 2.380E-01 
Iv= 14 1.674E+00 
iv= i3 3.908E+00 
im= 12 3.687E+00 
Ivy= 11 2.532E+00 
IY= 10 1.5389E+00 
IYy= 9 1.613E+00 
Iv= 8 1.759E+00 
IvY= 7 4.607E+00 
IY= 6 6.260E+00 


-2.960E-01 
-3.327E-01 
-7.082E-01 
-1.774E+00 
-2.398E+00 
-3 .835E+00 
-6 .009E+00 
-6.675E+00 
-7.184E+00 
-8.161E+00 
-8.362E+00 
-7.687E+00 
-7 .516E+00 
-6.517E+00 
=2 -522ueoO 
ee 


-9.456E-01 
-3.208E+00 
-5.263E+00 
-1.Q077E+0O1 
-1.197E+0O1l 
-1.063E+01 
-9.535E+00 
-9.226E+00 
-9.502E+00 
-6.387E+00 
-5.896E+00 
-5.322E+00 
-4.638E+00 
-2.466E+00 
eee 


lees OL 
Peli cle O nl 
-/Zoer00 
Sonne 
-637E+00 
-448E+00 
- 3468100 
»Z05En00 
. S608 + 00 
«92 5E00 
»714E+00 


MWIRrRRNWWERN OO 


MMW OOF WW ND 


.746E-01 
M7 5oE-0l 
~-461E-01 
-281E+00 
. 306E+00 
-416E+00 
-579E+00 
~-124E+00 
~392E+00 
. 380E+00 
-2606E+00 
-L67E+00 
-805E+00 
.056E+00 
-635E+00 
~313E+00 


-O22E+00 
.876E+00 
-456E+00 
-919E+00 
-I89E+00 
-882E+00 
~-341E+00 
~732E+00 
.770E+00 
-563E+00 
~-136E+00 
-990E-01 
.987E-02 
~167E+00 
2250E=01 


~739E-01 
.973E-01 
.582E+00 
~293E+00 
.043E+00 
.891E+00 
~597E-O1 
-O85E-01 
5Z1E=-01 
.758E+00 
~O26E+00 
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-6.328E-01 
-6.196E-01 
-1.135E+00 
-2.570E+00 
-4 .304E+00 
-6.973E+00 
-8.903E+00 
-9.284E+00 
-9.235E+00 
-1.074E+01 
-9.280E+00 
-8.019E+00 
-7 .572E+00 
-7.153E+00 
-6.075E+00 
Pe 


-9.759E-01 
-2.413E+00 
-3.610E+00 
-6.917E+00 
-7.933E+00 
-6.934E+00 
-5.091E+00 
-4.318E+00 
-4.646E+00 
1.539E-01 
2.705E-Ol1 
Z2.207E-01 
7.894E-01 
1.455E+00 
4s 


.334E-01 
.996E-01 
~-356E+00 
.601E+00 
.351E+00 
~-420E+00 
-779E-01 
-026E+00 
~-357E+00 
.665E+00 
»L33E+00 


UWNRrOr’rNNEFRNDND 


-7., 
~712E+00 
-447E+00 
~032E+00 
~325E+00 
-749E+00 
.602E+00 
-318E+00 
~-296E+00 
.964E-01 
.309E-02 
~LS7E-01 
-288E-01 
.781E-01 
~314E-01 


= 
2 
-4 


PWNYORPRRRPRONE 


-566E-01 
-716E-01 
.184E+00 
.555E+00 
-183E+00 
.956E+00 
.628E+00 
-O98E+00 
.288E+00 
ssole+OO 
.369E+00 
.688E+00 
~341E+00 
SOOLETOO 
-227E+00 
eee 


OSE =O) 


-.849E-01 
.834E-01 
-852E-01 
-681E+00 
.561E+00 
.466E+00 
-684E+00 
.984E+00 
7 nO 0 
.084E+00 
»-2L6E+00 


Be 5 1.186E+01 
ITy= 4 1.692E+01 
in— 3 2.589E+01 
IY= 2 2.428E+01 
hs 1 2-421E+01l 
IX= 1 
IY= 16 1.389E-01 
Ty= 15 2.289E-01 
IY= 14 7.295E-01 
fy 13 1.224E+00 
ITY= 12 1.204E+00 
iy= 11 1.259E+00 
Iv= 10 1.484E+00 
IY= 9 1.681E+00 
IY= 8 1.83ZE+00 
IY= 7 2.-336E+00 
IY= 6 2-85ZE+t00 
i= 5 3.34ZE+00 
IYy= 4 3.680E+00 
IyY= 3 3.901E+00 
IY= 2 4.070E+00 
IY= 1 4.077E+00 
TX= 6 
FIELD VALUES OF KE 
IvY= 16 3.551E+00 
IY= 15 3.571E+00 
IY= 14 1.800E+00 
IY= 13 1.336E+00 
IY= 12 2.064E+00 
ie 11 3.197E+00 
IY= 10 4.266E+00 
IY= 9 4.931E+00 
IYy= 8 5.536E+00 
IY= 7 1.260E+01 
IY= 6 1.486E+01 
i= 5 129975+01 
IY= 4 ZeUrsn+0 1 
IyY= 3 2.417E+01 
i— 2 2.841E+01 
Iv 1 Z2.950E+00 
IX= 1 
IY= 16 4.756E+00 
iea= 15 3. /00E+00 
IY= 14 &4.652E+00 
ITY= 13 6.297E+00 
T= 12 7.359E+00 
IY= 11 8.988E+00 
IY= 10 9.925E+00 
Iv= 9 9.921E+00 
t= 8 9.916E+00 
i= 7 1.080E+01 
IY= 6 150205701 
i= 5 9.485E+00 
i= 4 9.297E+00 
i— 3 9.416E+00 
a= 2 9.874E+00 
Iy= 1 1.005E+01 
IX= 
FIELD VALUES OF EP 
m= 16 Gee 23E700 
P= 15 Z 206 Lt+O0 
Iy= 14 9.861E-01 
rr 3 8.919E-01 
IY= 12 2.312E+00 
Iv= 1l 4.499E+00 
IY= 10 6 .645E+00 
Iy= 9 8.096E+00 
Iv= 8 9.374E+00 
Iy= 7 2.664E+01 
ry= 6 3.244E+01 
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NP RPP RRP RANMNWNrR HE WW 
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-215E+01 
-073E+01 
-068E+01 
-220E+01 
-280E+01 


-547E+00 
-248E+00 
-718E+00 
Po SE +00 
.199E+00 
-402E+00 
»322ZE+00 
-456E+00 
~335E+00 
-404E+01 
roo 2ZE+01 
-666E+01 
415E+01 
~-414E+01 
.691E+01 
-637E+00 


.7/Z1E+00 
-942E+00 
.007E+00 
-275E+00 
-441E+00 
.990E+00 
~542E+00 
-237E+01 
~439E+01 
.123E+01 
-666E+01 
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Ww 
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~l141E+01 
~-428E+01 
-590E+01 
-/37E+01 
./88E+01 


no 29E+00 
-334E+00 
-579E+00 
,oo2ZeE+O00 
-206E+00 
7655700 
-570E+00 
.O19E+01 
POsZE+O1 
7o57E+01 
-248E+01 
.O0S5E+01 
-/35E+00 
nOS7E+00 
- L40E+01 
.705E+00 


.640E+00 
-653E+00 
-279E+00 
.880E+00 
. 706E+00 
-LO6E+01 
.864E+01 
.308E+01 
-336E+01 
~185E+01 
. 765E+01 


on 


OY COOORPRRRrrORPWNE Ww 


NWNNNE WOON F OV 


PRR 
if 


ses +00 
OS) Ie ON 
lev E+O] 
sZt E+O1 
~-249E+01 


-6Z21E+00 
-950E+00 
-038E+00 
-625E+00 
.996E+00 
-326Et+00 
2O55E+01 
~-184E+01 
- LO6E+01 
-246E+01 
.0O49E+01 
.854E+00 
-141E+00 
Hio7E+O0 
063E+01 
.S48E-01 


-7 2368700 
-S14E+00 
-335E+00 
poz Ze +00 
-129E+00 
,910E+01 
-504E+01 
.889E+01 
fo? 2E +O | 
eee O 1 
ote +l 


OonomMN 
wn 


NNMNNNEF WUW ~) 


WOWOOMDMOWOrOWOUOOAMNWHN Ww 
in 


eeoen +00 
-874E+00 
wo on +00 
zon *O0 
oe +00 


.996E+00 
-446E+00 
.286E+00 
p20USE+00 
-431E+00 
-655E+00 
.943E+00 
.917E+00 
-697E+00 
.O90E+01 
. 760E+00 
.027E+00 
.868E+00 
PieizbsO0 
.867E+00 
po25L700 


-445E+00 
-408E+00 
-Z261LE+00 
-941E+00 
Poole tO 1 
.049E+01 
~407E+01 
PIOUETO 1 
.298E+01 
-680E+01 
Poz/et+O)] 


B= 15 
Ty= 4 
i= & 
n= 2 
rie & 
TX= 
ive 16 
iS 
IY= 14 
b= 13 
TY= 12 
TY= ll 
Ey= 10 
Iy= 9 
Y= 8 
m= 7 
IY= 6 
i= 6 
Iy= 4 
IY= 3 
m= 2 
Iy= 1 


hm A 


ee ey te ee eee 
ao 


eo oLETOL 
oon TUL 
-432E+01 
noo. EOL 
ssoeetO1 


-564E+00 
.995E+00 
-410E+00 
-387E+01 
7 29E+01 
-228E+01 
~505E+01 
-506E+01 
-509E+01 
-SL6E+01 
-629E+01 
-414E+01 
~358Et+01 
-379E+01 
~482E+01 
~539E+01 


FIELD VALUES OF Hl 
IY= 16 ZeOZELUS 


TY= 15 
IY= 14 
IY= 13 
TY= 12 
IY= 11 
[Y= 10 
m= 9 
IY= 8 
Iy= 7 
IY= 6 
IY= 5 
Iv= 4 
i= 13 

2 

uh 


IY= .962E+05 
IYy= ~962E+05 
IX= 

IY= 16 -962E+05 
IY= 15 .964E+05 
IvY= 14 ~962E+05 
f= 13 -962E+05 
IY= 12 ~962E+05 
IY= ll -962E+05 
IvY= 10 ~962E+05 
Ivy= 9 ~962E+05 
Ivy= 8 -962E+05 
IYy= 7 2624705 
IY= 6 ~-962E+05 
i= 5 .962E+05 
Ivy= 4 ~962E+05 
Iy= 3 ~-962E+05 
IY= 2 -962E+05 
a i -962E+05 
FIELD VALUES OF TMPIl 
IY= 16 2.951E+02 
= Ss 3.040E+02 
ITY= 14 2.994E+02 
i= 13 2.969E+02 
Iv= 12 2.966E+02 
eis) bi 2.964E+02 
Ivy= 10 2.961LE+02 
IY= 9 2.961E+02 
IYy= 8 2.960E+02 
IY= 7 2.956E+02 
IY= 6 2.956E+02 


MND PM MDID IPP PO POPE PO PO NN NNMNNNMNNNNNNN NW WY 
t 


~O52E+05 
O06E+05 
-980E+05 
~978E+05 
~97T6E+05 
~973E+05 
973E+05 
~972E+05 
.968E+05 
.968E+05 
.966E+05 
~965E+05 
.963E+05 


Nh Wh PhO LW 
N 


NON HO HONDO NNN OPO 


2597 barO | 
.781E+01 
~593E+01 
~L7OE+0O1 
~156E+01 


~.962E+05 
.042E+05 
-997E+05 
974E+05 
97TZ2E+05 
~970E+05 
~967E+05 
~966E+05 
.966E+05 
.964E+05 
~964E+05 
-963E+05 
~962E+05 
~-962E+05 
~962E+05 
.962E+05 


-950E+02 
FOZ 2ErO2 
“Jose+o2 
~962E+02 
-960E+02 
.958E+02 
.955E+02 
.955E+02 
~-954E+02 
~952E+02 
~-952E+02 


SiN eS 
WW 


NO BON PO BO HON ND DO FO 


.990E+01 
-438E+01 
-440E+01 
pLOGe+Ol 
Pzie rol 


7 Jo2E+05 


~951E+02 
~995E+02 
.965E+02 
~g5ge+02 
s2o26t02 
~951E+02 
.950E+02 
~950E+02 
.950E+02 
~950E+02 
~950E+02 
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1.881E+01 
lige ol 
1.869E+01 
2.205E+01 
opie + OR 


2.962E+05 
2.982E+05 
2.967E+05 
2.963E+05 
2.962E+05 
2.962E+05 
2.962E+05 
2.962E+05 
2.962E+05 
2.962E+05 
2.962E+05 
2.962E+05 
2.962E+05 
2.962E+05 
2.962E+05 
4 geeah oe 


~951E+02 
oo OE +O 2 
-955E+02 
oo e +0 2 
-J50EF02 
-950E+02 
-950E+02 
.950E+02 
-950E+02 
-950E+02 
~.950E+02 
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Zeon +o 1 
2 eet OL 
2.162E+01 
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eee ot 


2.962E+05 
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~.950E+02 
.955E+02 
-951E+02 
~-950E+02 
-950E+02 
~950E+02 
-950E+02 
-950E+02 
-950E+02 
-950E+02 
-950E+02 
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ae) Z2.954E 402 
IY= 4 Zee nO02 
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m= 2 2.950B+02 
I¥= | 2.-950B+402 
ix= il 
IY= 16 Zeosun+OZ 
m= 15 2-928 +02 
IY= 14 Zen +OZ 
IY= 13 Pee oUE +O 2 
IY= 12 2. ISOE +02 
E- 11 2e9 502102 
IY= 10 2.950E7 O02 
Iy= 9 2a s50EO02 
IY= 8 Z269e0e+02 
a 7 2e9 5002 
IY= 6 2. o50R +02 
B-— 5 Zea Ont OZ 
Iv= 4 A. GOETZ 
iy= 3 Zoo @E +O 2 
TY¥= 2 2. 950E4+02 
Ivy= 1 2.95002 
TA= 6 
PleLD VALUES OF" RHO] 
IY= 16 I. Zee+00 
TY= 15 Li tesE+00 
IY= 14 Lerose +00 
i= 13 Lose +00 
IY= 12 1 tee t00 
m— Ii 1. Peer eo0 
Py= 10 1. 19427-00 
Iy= 9 1. 295E400 
m—- 8&8 Lee +00 
IYy= 7 Te +00 
Iy= 6 Tete ErFOO 
m- 5 Pe 700 
IY= 4 lL eioer +00 
t7= 3 1 .SSE+00 
re 2 iden r00 
im 1 ie LO 9E+00 
ITA= 1 
Ty= 16 Peau lLe +00 
r— 15 LaZ20E +00 
TY= 14 T=Zoret00 
bee 3 ee Oe OO 
Ee I2 200 00 
ee 1 1. 200500 
IY= 10 1.200800 
Iv= 9 1 2Z00E+00 
IY= 8 Le200E+00 
=— 7 LeZzouE +00 
i—- 6 eZee +00 
m— 5 Ie 2008 +00 
Iy= 4 Te 2v0ETOO 
IY= 3 ee OE tO 0 
eo- 2 I2200E +00 
Iy= 1 Pee OOETO0 
TA= 6 
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~951E+02 
~951E+02 
.950E+02 
-950E+02 
-950E+02 


»201E+00 
- L69E+00 
-187E+00 
~-195E+00 
~1L95E+00 
-LS96E+00 
-197E+00 
~1L97E+00 
-197E+00 
~-L98E+00 
. LO98E+00 
- L98E+00 
.198E+00 
. LO8E+00 
-1L99E+00 
.199E+00 


~950E+02 
~-950E+02 
-950E+02 
~950E+02 
-950E+02 
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.183E+00 
~195E+00 
-199E+00 
-199E+00 
.199E+00 
-199E+00 
-199E+00 
-L99E+00 
.L99E+00 
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IY= 11 3.335E+02 
IY= 10 3.442E+02 
Iy= 9 3.535E+02 
Ivy= 8 3.594E+02 
IY= 7 3.799E+02 
IY= 6 4.239E+02 
Iy= 5 4.667E+02 
Ivy= 4 5 .939E+02 
Iy= 3 8.104E+02 
IvY= 2 6.146E+02 
Ivy= 1 6.043E+02 
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Iv= 16 3.453E+02 
Tv= 15 3.426E+02 
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-492E+02 
~382E+02 
~304E+02 
~349E+02 
~367E+02 
~419E+02 
~503E+02 
.555E+02 
~736E+02 
WO76E +02 
-377E+02 
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.479E+02 
.890E+02 
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3.499E+05 
3.490E+05 
3.366E+05 
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3.303E+05 
3.315E+05 
3.349E+05 
3.402E+05 
3.436E+05 
3.554E+05 
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~-476E+02 
.352E+02 
~-205E+02 
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~-476E+05 
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~246E+05 
-259E+05 
~-266E+05 
~289E+05 
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~349E+05 
~-430E+05 
~559E+05 
~667E+05 
-803E+05 
-834E+05 
~630E+05 
-487E+05 


ps See ees 
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~462E+02 
~325E+02 
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IY= 13 
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~-263E+02 
~L67E+02 
.166E+02 
~-168E+02 
. LB0E+02 
-L99E+02 
21L1LE+02 
~-253E+02 
~307E+02 
.349E+02 
~373E+02 
~-349E+02 
~271E+02 
»261E+02 
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-068E+00 
-052E+00 
.045E+00 
-028E+00 
~OO01E+00 
.843E-01 
~308E-01 
~347E-01 
-590E-01 
.073E-01 
~OO2E-01 
~/96E-01 
.939E-01 


-O26E+00 
.034E+00 
~-O86E+00 
»118E+00 
-L118zE+00 
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-113E+00 
- LO6E+00 
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.O70E+00 
~056E+00 
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.084E+00 
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.013E+00 
-014E+00 
.047E+00 
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~057E+00 
-051E+00 
~035E+00 
~O1LOE+00 
.947E-01 
-463E-01 
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-O70E-01 
~-134E-01 
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-284E-01 
-614E-01 
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.017E+00 
.019E+00 
.O056E+00 
.O084E+00 
.O76E+00 
.072E+00 
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~044E+00 
~.033E+00 
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-Ol1E-01 
~417E-01 
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~.772E-01 
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-023E+00 
.065E+00 
.095E+00 
-O90E+00 
.O088E+00 
.O80E+00 
.068E+00 
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.035E+00 
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~662E-01 
~310E-01 
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-O18E+00 
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-O22E+00 
-O28E+00 
.073E+00 
~ LO5SE+00 
»LOZE+O0 
. LOOE+00 
.094E+00 
.084E+00 
.078E+00 
~-058E+00 
-O28E+00 
~O05E+00 
~8l19E-01 
.814E-01 
~024E+00 
.042E+00 
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~823E+01 
~987E+01 
.924E+01 
-675E+00 
56 tse-01 
.O51E+01 
~.683E+02 
~C96E+02 
~771LE+O2 
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550 ZSLAB NO= 13 ITERN NO= 
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~-615E+t00 -8.031E+00 -7.528E+00 
~636E+01l -1.534E+01 -1.441E+01 
~612E+O1 -3.407E+O1l -3.239E+01 
.468E+01 -5.260E+01 -5.119E+01 
~LOZE+O1L -5.308E+01 -5.520E+01 
~322E+t01l -5.437E+01l -5.655E+01 
434E+O01 -4.802E+01 -5.308E+01 
-2l16E+O1l -3.185E+01 -4.241E+01 
-458E+01 -2.525E+01 -3.787E+01 
241Et+01 1.848E+0O1 -7.679E+00 
-2OLE+02 1.397E+02 7.413E+01 
-OLLE+02 2.499E+02 1.417E+02 
~ 740E+02 5.141E+02 2.898E+02 


.056E+00 
.341E+01 
~045E+01 
~920E+01 
.589E+01 
. /63E+01 
~654E+01 
.O38E+01 
~798E+01 
~917E+01 
.983E+01 
.653E+01 
~-203E+02 


IY= 3 1.558E+03 
TY= 2 -2.975E+01 
Iy= 1 8.651E+02 
IX= i; 
IY= 16 -5.601E+00 
ITY= 15 -9.771E+00 
TY= 14 -2.164E+01 
Py= 13 <3225122701 
by= 12 =326532201 
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IY= 10 -3.741E+01 
IY= 9 -£3.626E+01 
ITY= 8 -3.574E+01 
IY= 7 -3.155E+01 
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TY= 11 -3.050E-01 
IY= 10 -2.917E-01 
IY= 9 -2.698E-01 
IY= 8 -2.546E-01 
IY= 7 2-1.797E-01 
IY= 6 3.662E-01 
IY= 5 1.049E+00 
Iv= 4 5 .550E+00 
Y= 73 1.914E+01 
IY= 2 5.054E+00 
TY=s9 1. 3.463E+00 
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TY= 15 -1.135E+00 
TY= 14 -2.507E+00 
ITY= 13 -3.709E+00 
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IY= 11 -4.785E+00 
IY= 10 -4.765E+00 
IY= 9 °4.682E+00 
IY= 8 2*4.704E+00 
IY= 7 -4.660E+00 
IY= 6 -4.089E+00 
ITY= 5 -+4.054E+00 
IY= 4 -2.030E+00 
IY= 3 4.330E+00 
IY= 2 21.100E+01 
IY= 1 -2.810E-01 
IX= 1 
IY= 15 -8.082E-01 
IY= 14 -1.732E+00 
IY= 13 -2.360E+00 
IY= 12 -3.061E+00 
IY= 11 -3.129E+00 
IY= 10 -3.213E+00 
IY= 9 -3.335E+00 
IY= 8 =32.3808700 
IY= 7 -3.412E+00 
IY= 6 -3.322E+00 
IY= 5 -3.299E+00 
Iv= 4 -2.879E+00 
IY= 3 -2.508E+00 
IY= 2 <-1.849E+00 
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-4.340E-01 
-9.398E-01 
-1.218E+00 
-1.151E+00 
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-1.057E+00 
-9.226E-01 
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-661E+00 
.593E+00 
-450E+00 
~349E+00 
~J/14E-01 
-865E-02 
.O41E+00 
.564E+00 
-578E+00 
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-299E+00 
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.373E+00 
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-2.264E+00 
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-541E+00 
~570E+00 
.560E+00 
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.047E+00 
-ZL0E+00 
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.833E+00 
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- 798E+00 
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.780E+01 
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~358E+01 
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-445E+01 
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~587E+02 
-420E+02 
.803E+02 
~283E+02 
-2/79E+03 
.814E+03 
-161E+03 
~374E+03 
.OL8E+03 
~-239E+03 
-226E+04 
~-423E+04 
.955E+03 
.885E+01 
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Pook oo 
-483E+00 
~925E-01 
-543E+00 
-641E+00 
-512E+00 
~-274E+00 
~L57E+01 
-298E+01 
~772E+01 
~7/45E+01 
.564E+01 
-248E+01 
-623E+01 
-589E+01 
~532E+01 
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-287E+00 
-450E+00 
.862E-01 
~-491E+00 
~533E+00 
~372E+00 
-025E+00 
~LOVE+0O1 
-230E+01 
.636E+01 
~-375E+01 
.989E+01 
.O66E+01 
-856E+01 
./96E+01 
-615E+01 


.873E+01 
~-888E+01 
-518E+02 
~174E+02 
.892E+02 
~112E+02 
.092E+03 
~501E+03 
.762E+03 
.661E+03 
440E+03 
~914E+03 
-179E+03 
.987E+03 
.790E+03 
.896E+00 
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~234E+00 
-414E+00 
.603E-01 
Fo 7UE+OO0 
~312E+00 
~L1S5E+00 
-683E+00 
.0O58E+01 
~L7LE+O1l 
-538E+01 
Pl63E+O1 
.678E+01 
497E+01 
.OS8E+01 
.978E+01 
.829E+01 
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-673E+01 
-447E+02 
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ey = lez 4.282E+03 
bye" ii 5.716E+03 
IY= 10 9.936E+03 
Iy= 9 1.863E+04 
IY= 8 2.516E+04 
IY= 7 5.117E+04 
IY= 6 1.413E+05 
Iy= 5 2.325E+05 
IvY= 4 5.137E+05 
IYy= 3 8.560E+05 
Iv= 2 4.171E+05 
Iy= 1 3.094E+03 
IX= 1 
IY= 16 7.960E+01 
IY= 15 9.063E+01 
IY= 14 1.754E+02 
IY= 13 5.392E+02 
EY= "12 1.7 2¢5403 
Iv= 11 2.159E+03 
IY= 10 3.261E+03 
Iy= 9 5 .O77E+03 
iy= 8 6.299E+03 
IY= 7 1.078E+04 
IY= 6 2.055E+04 
IY= 5 2.911E+04 
IYy= 4 4.161E+04 
IY= 3 4.591E+04 
IY= 2 3.539E+04 
iy= 1 3.338E+04 
IX= 6 
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IY= 15 3.557E+05 
IY= 14 3.454E+05 
TY= 13 3.397E+05 
IY= 12 3.468E+05 
ey = 3.494E+05 
IY= 10 3.560E+05 
IY= 9 32 650E +505 
IYy= 8 3.722E+05 
IYy= 7 3.928E+05 
IY= 6 4.334E+05 
Iy= 5 4.673E+05 
IYy= 4 5 .427E+05 
IY= 3 6.691E+05 
ryv= 2 5.885E+05 
IY= 1 5.900E+05 
IX= ii 
IY= 16 3.510E+05 
ITY= 15 3.488E+05 
IY= 14 3.332E+05 
my = 213 3.251E+05 
Iv= 12 3.266E+05 
iey= )1 3.273E2+05 
IY= 10 3.293E+05 
Iy= 9 3.35232+05 
Iy= 8 3.341E+05 
lY= 7 3.403E+05 
Iy= 6 3.490E+05 
Ty= 5 3.558E+05 
lTY= 4 3.620E+05 
IY= 3 3.626E+05 
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~-484E+05 
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~851E+02 
,O085b703 
.016E+03 
.603E+03 
~1L47E+04 
~495E+04 
-832E+04 
-624E+04 
»-Q1OE+05 
~779E+05 
e238E+05 
~-366E+05 
~1Z1E+02 


-538E+05 
.538E+05 
-418E+05 
~348E+05 
~392E+05 
-409E+05 
~452E+05 
~517E+05 
~-557E+05 
-690E+05 
~-919E+05 
-LOPE+05 
~383E+05 
-545E+05 
~-434E+05 
-381lE+05 
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i 


~174E+01 
.570E+01 
-O28E+02 
.961E+02 
.509E+03 
~214E+03 
. L104E+03 
~476E+03 
.O83E+04 
-975E+04 
e271E+04 
~338E+04 
.O28E+05 


- 200m O02 
& 


~530E+05 
-524E+05 
NS 2E+O05 
~316E+05 
- 346E+05 
~358E+05 
~390E+05 
~438E+05 
-468E+05 
-567E+05 
- 728E+05 
-855E+05 
~O2Z2E+05 
eO095E+05 
.990E+05 
~9O0GE+05 


8.087E+01 
9.367E+01 
1.916E+02 
6.276E+02 
2.-148E+03 
2.-723E+03 
4.234E+03 
6 .845E+03 
8.646E+03 


-420E+04 


~522E+05 
.509E+05 
~368E+05 
~289E+05 
~3LLE+05 
sSe2LETO5 
~347E+05 
-386E+05 
~-410E+05 
-491E+05 
-617E+05 
.716E+05 
-832E+05 
~870E+05 
~772E+05 
~721E+05 


J OJ WW WO) WW WO) LW WO) C0 
un 


aY= 1 


FIELD VALUES 0 
iz= 16 ac 


mY- 15 
iy= 14 
[Y= 13 
mi 12 
IY= 11 
m= 10 
m- 9 
Iy= 8 
a «67/ 
IY= 6 
Iy= 5 
Iy= 4 
Y= 3 
IY= 2 
i= 1 
ITA= 
IY= 16 
IY= 15 
IY= 14 
IY= 13 
IY= 12 
TYy= 11 
IY= 10 
IY= 
IY= 
IY= 
IY= 
IY= 
IY= 
IY= 
IY= 
[Y= 


Po OJ Pf U10) ~1 0 0 


IX= 
BEELD VAL 
I 


Y= 16 
IYy= 15 
[Y= 14 
Ex= 13 
IY= 12 


oe 
6 


WOUWDWUORP PRR RP RRP RRR DNHNUN~)TDWOWOWOWORPRPRPRPHP RPG WWWWWWWWWWWWWWWW WOOP BW WWW WWW Wl 
oo ° OV f= 


5525205 


38E+02 


-542E+02 
-440E+02 
~383E+02 
~-454E+02 
-450E+02 
-546E+02 
-646E+02 
- 708E+02 
~-913E+02 
~SLITEFTO2 
-655E+02 
-405E+02 
-664E+02 
~-862E+02 
~877E+02 


-496E+02 
-474E+02 
~31L9E+02 
238E+02 
~-253E+02 
-2600E+02 
~-279E+02 
» 309E+02 
~328E+02 
- 569E+02 
»-476E+02 
.544E+02 
-606E+02 
~6O11LE+02 
.548E+02 
Sooo TOZ 


Boer enn l 
00 


LE+00 


.O00E+00 
.029E+00 
.047E+00 
0258700 
.018E+00 
.989E-01 
~717E-01 
.556E-01 
.O059E-01 
-220E-01 
-625E-01 
.537E-01 
-398E-01 
.042E-01 
-O80E-01 


.O13E+00 
-O20E+00 
-067E+00 
.094E+00 
.O389E+00 
.086E+00 
.O80E+00 
-.070E+00 
.064E+00 
.045E+00 
.O19SE+00 
299 5h=0 1 
~825E-01 
~8llE-01 
.986E-01 


BT MPa 


10101 BB OD LO Ld Ld G) Ld WW Ld Ld 
No 


SAHAAIIDOWOOPRPRREEE 
NO 


.534E+02 
-539E+02 
~433E+02 
.373E+02 
-438E+02 
-462E+02 
~523E+02 
-614E+02 
~-670E+02 
-855E+02 
-L95E+02 
~-466E+02 
~-958E+02 
~310E+02 
~L62E+02 
-L16E+02 


.002E+00 
-O001E+00 
s0sZ25700 
.OS0E+00 
-030E+00 
-023E+00 
.005E+00 
-801E-01 
~-651E-01 
~-193E-01 
-464E-01 
.964E-01 
-214E-01 
-755E-01 
.956E-01 
~O27E-01 


BB BB Lo WW WW Wd WW WW Wo 


DWOO~ATNDWDMWOWOWORPRPRRRRRH 
J 


~-524E+02 
-524E+02 
~-405E+02 
~335E+02 
.379E+02 
.395E+02 
~-438E+02 
-503E+02 
-543E+02 
~.675E+02 
.904E+02 
.O85E+02 
-365E+02 
~527Et+02 
-416E+02 
. 363E+02 


-005E+00 
-005E+00 
.040E+00 
.062E+00 
-048E+00 
.043E+00 
.030E+00 
Polite +00 
.998E-01 
~.642E-01 
~O88E-01 
-694E-01 
~1L57E-01 
.884E-01 
.0O75E-01 
-176E-01 
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OJ) BoB OG) OW) 9 Ld Od LW Wd Ld Ws Ld OG 


WODWDOWWDOUWORPRPRPRPRPRPRREH 
fs 


-516E+02 
- 509E+02 
»379E+02 
. 303E+02 
.333E+02 
-345E+02 
~377E+02 
-424E+02 
454E+02 
Pa oSE+O 2 
So SEFOZ 
.839E+02 
.O06E+02 
.O79E+02 
-974E+02 
.890E+02 


-0O0SE+00 
.009E+00 
.048E+00 
~-072E+00 
.062E+00 
.059E+00 
.049E+00 
.034E+00 
~.O25E+00 
~971E-01 
.548E-01 
-240E-01 
~.870E-01 
~720E-01 
.947E-01 
~141E-01 


ld 0 Ld LY Ld Wd Ld Wo Wd WW WW Ws Ws Wd 
un 


WOWWWWDRPRPRPERHRPRPRHEHE 
un 


.508E+02 
-495E+02 
.354E+02 
-276E+02 
POSE 
. 308E+02 
~333E+02 
~372E+02 
ey sol 
~-478E+02 
~603E+02 
~701LE+02 
-S17E+02 
ODD E+O2 
»f57E+02 
- /06E+02 


NU OE+O0 
~O13E+00 
.056E+00 
-O0O81E+00 
.074E+00 
-O71ETOO 
.062E+00 
-050E+00 
.043E+00 
.018E+00 
.836E-01 
~577E-01 
.293E-01 
-2Z07E-01 
~445E-01 
.585E-01 


TY= 
IXx= 


i 


1). 002E7,00 
6 


KAAERKKKKKRKKKKRKKKRKRKKRKKRKKRKKKKKEKKEKKKRRKKKKRKRKRKKRRKKRKRKRKKKKE 


TIME STP= 


FLOW FIELD AT ITHYD= 


FIELD VALUES OF Pl 


Ty= 16 2.236E+00 
IY= 15 -9.074E+00 
IY= 14 -3.507E+01 
IY= 13 -6.000E+01 
IY= 12 -8.413E+01 
TY= 11 =-8.255E+01 
IY= 10 -6.384E+01 
IY= 9 -2.785E+01 
IY= 8 -1.335E+00 
IvY= 7 1.551E+02 
IY= 6 5 .847E+02 
TY= 5 1.154E+03 
Iy= 4 2.593E+03 
IY= 3 6.513E+03 
IY= 2 1.181E+03 
Ivy= 1 1.118E+03 
IXx= il 
IY= 16 1.126E+00 
IY= 15 -6.579E+00 
IY= 14 -2.500E+01 
IvY= 13 -4.612E+01 
IY= 12 -6.926E+01 
IY= 11 -7.282E+01 
IY= 10 -7.849E+01 
ITY= 9 -8.425E+01 
IY= 8 -8.701E+01 
IY= 7 -9.108E+01 
IY= 6 -9.580E+01 
TyYv= 5 -+9.742E+01 
IY= 4 -1.015E+02 
IY= 3 -1.132E+02 
ITY= 2 -1.380E+02 
IY= 1 -1.056E+02 
IX= 6 
FIELD VALUES OF Ul 
IY= 16 -5.286E-02 
TY= 15 -9.333E-02 
TY= 14 -2.229E-01 
TY= 13 =2.8018-01 
IyY= 12 -2.407E-01 
[Y= 1] =232772-01 
TY= 10 -1.858E-01 
IY= 9 -1.059E-01 
IY= 8 -4.503E-02 
Ivy= 7 2.2608E-01 
IY= 6 1.355E+00 
i= 5 2.6/5E+00 
Iv= 4 8.192E+00 
IY= 3 1.936E+01 
Ivy= 2 7.762E+00 
Ivy= 1 6.268E+00 
IXx= 1 
FIELD VALUES OF V1 
ITY= 15 -6.343E-01 
IY= 14 -1.622E+00 
IY= 13 -2.650E+00 
IY= 12 -3.520E+00 
IY= 11 -3.527E+00 
TY= 10 -3.455E+00 
IY= 9 -3.242E+00 
IY= 8 -3.161E+00 


1 SWEEP NO= 


=a 
=3 
= 
= 
= 
= 
=5 
= 
= ll 


be 
2s 
Se 
g. 
i 
9 
oF 

Z 


1, 


OOOE+01 


~-677E+01 
~2UZET0Z 
-549E+02 
.863E+02 
.979E+03 
~243E+03 
-L126E+03 
Nee 


335E-01 
153E-01 
699E-01 
460E-01 
656E-01 
052E-01 
423E-01 
856E-01 
136E-01 
230E-01 
172E+00 
767E+00 
O46E+00 
LIV E+FOL 
415E+00 
774E+00 


~-302E-01 
-610E+00 
-61LLE+O0 
-471E+00 
-474E+00 
-401E+00 


183E+00 
O90E+00 


IZ= 14, 


2.057E+00 
-1.008E+01 
-3.876E+01 
-6 .605E+01 
-9.118E+01 
-8.974E+01 
-7.232E+01 
“4, 
-1 
1 
4 
8 
J 
Z 
1 
1 


550 ZSLAB NO= 


ISWEEP= 550, 


Z2eS5E+00 
-9.541E+00 
-3.,805E+01 
-6.967E+01 
-1.092E+02 
-1.128E+02 
-1.119E+02 
-1.072E+02 
-1.023E+02 
-5.161E+01 
.894E+01 
-L90E+02 
~561E+02 
- 789E+02 
-526E+02 
-LO9E+02 


H NO 0. UNO OV 
Ld 


-1.621E-01 
-4.174E-01 
-1.030E+00 
-1.317E+00 
-1.184E+00 
-1.133E+00 
-9.797E-01 
=7 ote =01 
-5.646E-01 
4.491E-02 
1.489E+00 
2.821E+00 
5 .880E+00 
7.989E+00 
7.056E+00 
7.391E+00 
3 
-5.558E-01 
-1.460E+00 
-2.382E+00 
-3.290E+00 
-3.313E+00 
-3.306E+00 


-3.243E+00 
-3.202E+00 
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14 ITERN NO= if 

ISTEP= 

2.246E+00 2.354E+900 
-9.047E+00 -8.500E+00 
-3.766E+01l -3.730E+01 
-7.348E+01 <-7.738E+01 
-1.234E+02 -1.356E+02 
-1.308E+02 -+1.460E+02 
-1.404E+02 -1.630E+02 
“1.51@E+O02 <-1.625E7102 
-1.559E+02 -1.930E+02 
-1.489E+02 -+2.102E+02 
-“1.250E+t02 <-2.307E+02 
-8.656E+01 -2.449E+02 
8.401E+00 -2.436E+02 
Z2.087E+O1l -2.728E+02 
-8.2/6E+01 “=se303Er02 
1.308E+0Z <2.352E+t02Z 

4 > 
-1.739E-01 <-1.649E-01 
-5.122E-01 -5.477E-01 
-1.298E+00 -1.436E+00 
-1.759E+00 -2.039E+00 
-1.794E+00 -2.259E+00 
-1.779E+00 -2.282E+00 
-1.690E+00 -2.255E+00 
-1.510E+00 -2.146E+00 
-1.383E+t00 -2.060E+00 
-9.052E-01 <-1.683E+00 
1.864E-01 -8.099E-01 
1.171E+00 -1.582E-02 
3. 2352700 1.492E+00 
4.363E+00 2.323E+00 
3.942E+00 1.782E+00 
4.089E+00 1.834E+00 

4 5 
-4.884E-01 -4.287E-01 
-1.324E+00 <-1.197E+00 
-2.165E+00 <-1.947E+00 
-3.072E+00 -2.795E+00 
-3.109E+00 -2.843E+00 
-3.147E+00 -2.912E+00 
-3.190E+00 <-3.016E+00 
-3.190E+00 -3.045E+00 


TS ae 


IY= 7 -+2.997E+00 
IY= 6 -1.484E+00 
IY= 5 -1.226E+00 
Iy= 4 1.370E+00 
te 3 4.215E+00 
IY= 2 -1.189E+01 
ITY= 1 -3.387E+00 
IX= I 
TY= 15 -3.462E-01 
IY= 14 -1.0135&+00 
TY¥= 13 -1.612E+00 
ITY= 12 -2.299E+00 
TY= 11 -2.375E+00 
IY= 10 -2.476E+00 
IY= 9 -2.629E+00 
IY= 8 -2.694E+00 
IY= 7 -2.750E+00 
IY= 6 -2.731E+00 
ITY= 5 +2.760E+00 
IY= 4 -2.487E+00 
IY= 3 -2.299E+00 
IY= 2 <-1.829E+00 
a 1 i ae 
FIELD VALUES OF Wl 
IY= 16 -2.707E+00 
ITY= 15 -1.563E+00 
Iv= 14 9.543E-01 
f= 13 3.767E+00 
y= 12 7.218E+00 
IY= ll 8.162E+00 
IY= 10 1.022E+01 
IYy= 9 1.316E+01 
Iy= 8 1.504E+01 
iy—= 7 2. 147E+01 
IY= 6 3.734E+01 
ITy= 5 5.125E+01 
Ty= 4 8.627E+01 
Iy= 3 125308402 
IY= 2 1.279E+02 
Ivy= 1 ISS 7Zet02 
IX= i 
ITY= 16 -2.468E+00 
TY= 15 -1.434E+00 
IY= 14 1 .029E+00 
fy= 13 3.793E+00 
IY= 12 7.195E+00 
IY= 11 8.054E+00 
TY= 10 9.728E+00 
Iy= 9 1.169E+01 
IvY= 8 1.283E+01 
IY= 7 1.644E+01 
IY= 6 2.197E+01 
Iy= 5 2.646E+01 
IY= 4 3.264E+01 
ia 3 3.668E+01 
im 2 3.724E+01 
ha 1 ao 


PLELD VALUES OF KE 


IY= 
Iv= 
IY= 
Iy= 
IY= 
IY= 
IY= 
iY= 
IY= 
Iv= 


IL 


OWWNRPReR PNF 


~.094E+02 
-230E+02 
-L37E+02 
1 2UE Og 
-084E+03 
~287E+03 
pi eze+0 3 
-594E+03 
.139E+03 
eet e+O3 


PONDER BNR EH 


MPRiO~A) BWR FH O~I~AWWFRND 
NO 


.909E+00 
. 306E+00 
~978E-01 
. 200E+00 
2555-01 
-625E+00 
Seay 


.719E+00 
-594E+00 
.O96E-01 
-686E+00 
.052E+00 
.966E+00 
-926E+00 
-268E+01 
-442E+01 
.O30E+01 
.396E+01 
»D04E+01 
.2/2E+01 
.695E+01 
.002E+02 
.044E+02 


1092002 
~-226E+02 
Sez e+OZ 
.647E+02 
.046E+03 
-236E+03 
~695E+03 
.426E+03 
.909E+03 
-603E+03 


WNRRRO PN Fr 


' 3 
DADAM WNOR RF WO ~)~)0I WF fh 
iu 


PEOZE+00 
» LS6E+00 
-996E+00 
-842E-01 
-430E+00 
-688E+00 
ene 


Poo 7E200 
.563E+00 
-210E-01 
moore +00 
.023E+00 
~910E+00 
e700 
-214E+01 
esere7+O | 
.846E+01 
~837E+01 
.666E+01 
232Z5E+01 
ogg e FOL 
.874E+01 
.982E+01 


mOeSe +02 
~209E+02 
~049E+02 
.333E+02 
.145E+02 
.067E+03 
.421E+03 
.954E+03 
.299E+03 
-495E+03 
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i 4 
POOP WR RRWO~IDWWFr ho 


NRE OOPRRH 


~147E+00 
-638E+00 
-565E+00 
-825E+00 
-014E+00 
-568E+00 
-140E-01 


-591E+00 
~530E+00 
.284E-01 
.674E+00 
.978E+00 
.842E+00 
»578E+00 
175E+01 
.304E+01 
~7Z27E+01 
nego +O 1) 
129E+01 
-231E+01 
~O5S5E+01 
. LO3E+01 
.888E+01 


~O75E+02 
~193E+02 
~982E+02 
-087E+02 
~-239E+02 
.516E+02 
Zoe Fos 
.649E+03 
-910E+03 
PoesEtTUS 


NRE HO~AWRHe 


DPOB UB WN MRR RFHWO~IDWWrH fo 
wn 


.047E+00 
.785E+00 
. 764E+00 
-240E+00 
eZ roe+O0 
.183E+00 
eee ot 


~538E+00 
-494E+00 
~423E-01 
.660E+00 
~954E+00 
.805E+00 
-489E+00 
~154E+01 
.275E+01 
.664E+01 
~322E+01 
.861E+01 
~716E+01 
Poe +O] 
~32Z2E+01 
.209E+01 


.068E+02 
~175E+02 
Pousn+0 2 
~874E+02 
.602E+02 
~724E+02 
meZOE+O3 
.464E+03 
.679E+03 
-411E+03 


Iy= 6 1.020E+04 
Iy= 5 1.477E+04 
IY= 4 2.664E+04 
Iy= 3 4.607E+04 
IY= 2 2 WISET04 
Iy=- Tea oeror 
IX= A 
IY= 16 TOs Ee 
IY= 15 1. LEGER BOZ 
IY= 14 1.804E+02 
Y= 3 3 OSCE +O Z 
= eZ 6a OPE OZ 
Ey= 11 TeGa7buoZ 
ITY= 10 92 70SE 02 
r= @2 lL 2ooenoS 
Iy= 8 1.404E+03 
i= 7 1. Soon 0S 
IY= 6 Z2.670E803 
i> 5 JeOSoedO Ss 
IY= 4 4 .530E+03 
i= 2 4.760E+03 
IY= 2 3.942203 
re ft 3, S0DE+03 
IA= 
PIELD VALUES@@ES Ee 
IY= 16 lL. LOSE TO Z 
ie 15 ee eee nO 2 
IY= 14 3.438E+02 
TY= 13 T27oe 03 
IY= 12 479618703 
Ex ii 6.46 28+03 
IY= 10 1.078E+04 
IY= 9 1.957E+04 
IY= 8 2.613E+04 
IY= 7 5.192E+04 
-- 6 1.402E+05 
i- © ZeeoeEt05 
Iy= 4 &.767EF05 
i= 3 622562705 
IY= 2 3c enn +05 
IY= 1 See et03 
IA= 1 
IY= 16 Le U66E 702 
IY= 15 Zonet 2 
IY= 14 Zeon 
Iv= 13 6,0708702 
TY= i2 2.250E103 
BY =i Zeloounus 
IY= 10 Seozce sos 
Iz= 9 5.6252703 
IY= 8 7 Oe tEO S 
IY= 7 1.169E+04 
IY= 6 2.1/CE+04 
IY= 5 3.054E+04 
IY= 4 4.360E+04 
in- 38 4.818E+04 
i= 2 3.762ZE+04 
Iv= 1] 325758704 
IX= 6 
PIELD. VALUES GR Hi 
IY= 16 3.626E+05 
Iv= 15 3 OO9E EOS 
IY= 14 S7oeen 05 
IY= 13 S47 n0 o 
Iv= 12 3.546E+05 
Iy= 11 Se oe O 5 
Ty= 10 3.634E+05 
i= S37 2G FOS 
Iv= 8 Sageoe to 5 
IY= 7 Seo oe TO5 


.694E+03 
-216E+04 
.949E+04 
-406E+04 
.501E+04 
-462E+01 


Hb NO F F OO 
N 


1.106E+02 
1.346E+02 
3.399E+02 
1.233E+03 
4.676E+03 
6.033E+03 
9.885E+03 
1.750E+04 
2.304E+04 
4.444E+04 
1.113E+05 
1.727E+05 
3.235E+05 
4.217E+05 
2.576E+05 
‘ae 


~-623E+05 
~606E+05 
~504E+05 
-460E+05 
~531E+05 
-555E+05 
~6513E+05 
.699E+05 
»toage+O5 
~.924E+05 


Lo Wo Wo Wd Wd WY Co) Wd Wd 


ly Wo Wo Wo Wd WW WW 


N wrk OO Or 
ly 


PRR ODMR ~1B Wr Wr 
iy 


TEOOE +03 
-281E+03 
-221E+04 
-411E+04 
O18E+03 
~1L41E+01 


.093E+02 
-321E+02 
-228E+02 
.089E+03 
~720E+03 
~710E+03 
~402E+03 
~241E+04 
~597E+04 
~952E+04 
. 809E+04 
Ores OS 
~786E+05 
-199E+05 
~371E+05 
-989E+02 


~614E+05 
~593E+05 
-479E+05 
~425E+05 
~474E+05 
-491E+05 
Ooo E05 
-594E+05 
~632E+05 
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.789E+03 
~ LO3E+03 
Z81E+03 
-822E+03 
.609E+03 
»203E+03 
.794E+03 
~792E+03 
-685E+03 


pooeee= 


~574E+02 


ee 15 3.750E+02 
ivy= 14 1.090E+03 
fo 13 3.160E+03 
i= 12 8.341E+03 
IY= ll 9.679E+03 
IY= 10 1.302E+04 
Iy= 9 1.864E+04 
lYy= 8 2.217E+04 
m— 7 3.357E+04 
IyY= 6 6.485E+04 
m= 5 8.978E+04 
Iy= 4 1.532E+05 
Iy= 3 2.305E+05 
i= 2 1.174E+05 
IYy= 1 1.454E+03 
IX= 1 
IY= 16 2.564E+02 
Iivy= 15 3.543E+02 
IY= t4 9.480E+02 
IYy= 13 2.482E+03 
IY= 12 5 .602E+03 
IY= ll 6.358E+03 
i= 10 8 .OO00E+03 
Iy= 9 1: 025E+04 
Iy= 8 1.158E+04 
i= 7 1.576E+04 
IY= 6 2.292E+04 
Iy= 5 2.850E+04 
Iy= 4 3.495E+04 
i= 3 3.630E+04 
Iy= 2 3.098E+04 
IY= 1 2.992E+04 
IX= 6 
FIELD VALUES OF Hl 
IY= 16 3.915E+05 
Iy= 15 3.871E+05 
IY= 14 3.819E+05 
IY= 13 3.858E+05 
IvY= 12 3.964E+05 
IY= ll 3.987E+05 
Ivy= 10 4.034E+05 
Iy= 9 4.093E+05 
Iy= 8 4.123E+05 
Iy= 7 4.210E+05 
IY= 6 4.347E+05 
Iy= 5 4.446E+05 
Iy= 4 4.629E+05 
Iy= 3 4.953E+05 
ty— 2 5 .008E+05 
Iy= 1 5.051E+05 
IX= 1 
IY= 16 3.891E+05 
IY= 15 3.852E+05 
i= 14 3.769E+05 
i= 13 3.788E+05 
Iv= 12 3.860E+05 
Iv= 11 3.877E+05 
IY= 10 3.908E+05 
Ivy= 9 3.941E+05 
Iy= 8 3.959E+05 
Iv= 7 4.013E+05 
IY= 6 4.079E+05 
Ivy= 5 4.130E+05 
Iy= 4 4.182E+05 
Iy= 3 4.208E+05 
Iy= 2 4.200E+05 
Iv= 1 4.197E+05 
Ik= 6 
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NO 


BoP PoP BOB BBB WW WW WW 
NO 


~/34E+02 
.0O79E+03 
~ LOOE+03 
.075E+03 
~349E+03 
~-250E+04 
~770E+04 
.093E+04 
~L27E+04 
.801E+04 
~859E+04 
~-235E+05 
~394E+05 
-304E+04 
~432E+02 


~913E+05 
.868E+05 
~815E+05 
~852E+05 
~954E+05 
~977TE+05 
023E+05 
.O80E+05 
~L1LOE+05 
~194E+05 
~324E+05 
~417E+05 
~584E+05 
~757E+05 
.853E+05 
~-876E+05 


897E+02 


Hm 0200 00 OUP DFR FR O~INF WW 


or 


PP PB BBB BB lO dW WW to Wd 
lJ 


-660E+02 
.033E+03 
.859E+03 
-O60E+03 
~122E+03 
.065E+04 
-460E+04 
~ /02E+04 
~478E+04 
-246E+04 
~615E+04 
~.O12E+04 
.SO0E+04 
~IJ10E+04 
~-976E+02 


~904E+05 
~-862E+05 
-802E+05 
~830E+05 
~918E+05 
-938E+05 
~978E+05 
.O25E+05 
-O50E+05 
~124E+05 
~233E+05 
~314E+05 
~438E+05 
-538E+05 
-582E+05 
.591LE+05 


889E+02 


IES 


SoCo eos 


NOP ODOUR WORF WD ON WOW 
An 


pe Pees ee 


~609E+02 
.999E+02 
.696E+03 
. 388E+03 
ow +03 
~417E+03 
~-256E+04 
-446E+04 
.053E+04 
-295E+04 
~.263E+04 
.677E+04 
~O38E+04 
-O80E+04 
.665E+02 


~898E+05 
~858E+05 
~791LE+05 
~815E+05 
.893E+05 
~9L1LE+05 
.946E+05 
~986E+05 
~O09E+05 
~O75E+05 
~165E+05 
~234E+05 
~325E+05 
~387Et+05 
.393E+05 
~368E+05 


om 


NO WP BWR EF 00 OOM WO W 
in 


PPO BOB BB bo Wo WD) WW Wd Wd Wd Wd 
on) 


-576E+02 
77I2E+02 
~597E+03 
.O14E+03 
~854E+03 
~ 742E+03 
-1L46E+04 
. 309E+04 
~827E+04 
.S18E+04 
~593E+04 
~620E+04 
.841E+04 
~-567E+04 
~S18E+04 


~O95E+05 
~855E+05 
~782E+05 
-803E+05 
.878E+05 
~895E+05 
~928E+05 
.965E+05 
~985E+05 
.046E+05 
~125E+05 
.L86E+05 
~259E+05 
~302E+05 
.295E+05 
-277E+05 


879E+02 


Ty= 55 3.855E+02 
IY= 14 3.804E+02 
IY= 13 3.843E+02 
TY= 12 3.948E+02 
IY= 11 3.971E+02 
IY= 10 4.018E+02 
Iyv= 9 4.076E+02 
Ivyv= 8 4.107E+02 
iz= 4.193E+02 
IY= 6 4.330E+02 
k= 5 4.428E+02 
Iy= 4 4.610E+02 
Iy= 3 4.933E+02 
IY= 2 4.988E+02 
Iy= 1 5.031E+02 
IX= 1 
IY= 16 3.875E+02 
TY= 15 3.837E+02 
Iv= 14 3.754E+02 
Ivy= 13 3.773E+02 
IY= 12 3.845E+02 
IY= ll 3.861E+02 
IY= 10 3.892E+02 
IY= 9 3 ee 20 2 
IvY= 8 3.943E+02 
Iv= 7 3.997E+02 
IY= 6 4.063E+02 
IY= 5 4.113E+02 
Iy= 4 4.166E+02 
Ivy= 3 4.191E+02 
IY= 2 4.183E+02 
Ivy= 1 4.180E+02 
IX= 6 
FIELD VALUES OF RHO1 
Iv= 16 9.088E-01 
Ty= 15 9.190E-01 
IY= 14 9.313E-01 
Iv= 13 9.219E-01 
TY= 12 8.978E-01 
IY= ll 8.929E-Ol 
IY= 10 8.834E-01 
Iy= 9 8.724E-01 
IYy= 8 8.669E-01 
Iy= 7 8.531E-01 
IY= 6 8.330E-01l 
Iy= 5 8.210E-01 
Iy= 4 7.959E-01 
IYy= 3 7.278E-Ol 
IYy= 2 7.139E-01 
Iy= 1 7.083E-01 
IxX= 1 
IY= 16 9.144E-01 
mY= 15 9.234E-01 
IY= 14 9.437E-01 
TY= 13 9.388E-01 
IY= 12 9.214E-01 
IY= 11 9.175z£-01 
IY= 10 9.104E-01 
IY= 9 9.028E-0l 
Ivy= 8 8.987E-Ol 
f= 7 8.867E-01 
i= Fo 8.725E-01 
iy 5 8.620E-O01l 
Ivy= 4 8.512E-01 
i= 3 8.460E-01 
ine— 2 8.475E-O1 
IyYy= |] 8.481E-01l 
IAZ= 6 
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-Oeen=0 lL 
~468E-01 
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~L14E-01 
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3.842E+02 
3.7/6E+02 
3.799E+02 
3.878E+02 
3.895E+02 
3.930E+02 
3.971E+02 
3.993E+02 
4.058E+02 
4,.149E+02 
4.217E+02 
4.308E+02 
4.369E+02 
4.375E+02 
ose a 


9.1285-01 
9.220 -O1 
9esere-O1 
geoeoe—0] 
oe OL 
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JO le O 1 
azo eo 
8.881E-01 
8. 7405-01 
Been Ol 
8.445E-01 
S201 
ooo - Ol 
Sao OL 
= 


3.840E+02 
3.767E+02 
3.788E+02 
3.862E+02 
3.879E+02 
3.912E+02 
3.949E+02 
3.969E+02 
4.030E+02 
4,109E+02 
4.170E+02 
4.242E+02 
4.285E+02 
4.277E+02 
sper ee 2 


Ses 6k s01 
Sea 1 
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3 Ge Cee oe 
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8. 2535E-01 
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TOME STP= 1 SWEEP NO= 550 ZSLAB NO= 17 ITERN NO= 1 


FLOW FIELD AT ITHYD= 


iy a= dy LONEEF="550, ISTEP= i 
Peet) VALUES OF Pl 


IY= 16 4.164E+01 4.478E+01 4.398E+01 4.335E+01 4.277E+01 
TY= 15 1.860E+01 2.304E+01 2.278E+01 2.239E+01 2.195E+01 
TY= 14 -6.057E-01 6.623E+00 8.408E+00 9.779E+00 1.100E+01 
TY= 13 -1.357E+01 9.575E+00 1.287E+01 1.531E+01 1.762E+01 
TY= 12 -3.070E+01 2.455E+01 2.261E+01 2.134E+01 2.165E+01 
TY= 11 -3.314E+01 3.177E+01 2.725E+01 2.483E+01 2.429E+01 
TY= 10 -2.890E+01 5 .580E+01 4.266E+01 3.581E+01 3.236E+01 
IYy= 9 4.899E+01 1.534E+02 6.035E+01 4.518E+01 3.738E+01 
teow 8 5 .249E+01 1.697E+02 6.504E+01 4.527E+01 3.528E+01 
i 7 1.206E+02 2.750E+02 1.180E+02 7.706E+01 5.487E+01 
i= 6 1.963E+02 3.900E+02 1.719E+02 9.808E+01 6.190E+01 
Ty= 5 3.224E+02 5.426E+02 2-281E+02 1.258E+02 7.712E+01 
Iy= 4 4.290E+02 6.928E+02 2.508E+02 1.264E+02 8.230E+01 
Iy= 3 3.130E+02 4.059E+02 -1.354E+02 5.760E+01 6.294E+01 
TY= 2 3.683E+02 4.059E+02 1.383E+02 1.029E+02 8.407E+01 
iz— 1 4.029E+02 4.310E+02 lee S0E+0iZ 8.800E+01 6.505E+01 
IX= ll 2 3 4 5 

IY= 16 3.378E+01 

iz= 15 1.680E+01 

IY= 14 1.086E+01 

t= 13 2.083E+01 

f= 12 3.499E+01 

TY= ll 3.865E+01 

TYy= 10 4.588E+01 

Iy= 9 5 .264E+01 

Iy= 8 § ..535E+01 

lu 7 6.750E+01 

IY= 6 7.460E+01 

Tyv= 5 8.077E+01 

ly= 4 7.762E+01 

IY= 3 7.Q019E+01 

IY= 2 7.052E+01 

i ii 6.221E+01 

BeeeD VALUES OF Ul 

Iv= 16 4.045E-01 1.043E+00 1.399E+00 1.807E+00 2.017E+00 
ly 15 2.094E-01 7.736E-01 1.042E+00 1.312Z2E+00 1.435E+00 
TY= 14 1.950E-01l 7.116E-01 8.870E-0Ol 1.013E+00 1.059E+00 
IY= 13 4.788E-01 1.726E+00 2.065E+00 2.238E+00 2.278E+00 
ITY= 12 1.073E+00 3.781E+00 4.325E+00 4.455E+00 4.417E+00 
IY= 11 1.188E+00 4.222E+00 4.801E+00 4.922E+00 4.874E+00 
TY= 10 1.441E+00 5.211E+00 5.802E+00 5 .877E+00 5 .803E+00 
IYy= 9 1.814E+00 6.679E+00 7.075E+00 7.013E+00 6.873E+00 
i= 8 2.003E+00 7.381E+00 7.744E+00 7.619E+00 7.443E+00 
i” 7 2.509E+00 9.213E+00 9.592E+00 9.333E+00 9.088E+00 
IY= 6 3.119E+00 I Z28E+01 1.152E+01 1.096E+01 1.057E+01 
i= 5 3.506E+00 1.252E+01 1.267E+01 1.193E+01 1.148E+01 
ie 4 3.571E+00 1.237E+01 1.167E+0O1 1.100E+01 1.066E+01 
i= 3 1.706E+00 3.001E-04 5 .606E+00 7.976E+00 8 .464E+00 
IY= 2 2.480E+00 7.324E+00 7.836E+00 8.425E+00 8.521E+00 
i= 1 2.250E+00 2 GL aE LO S36 GEO BegeoEtOO Seeger 08 

S. i 

Paee LD VARUES OF V1 

TY= 15 5.579E+00 5 .586E+00 5 .564E+00 5 .566E+00 5.596E+00 
TY= 14 7.941E+00 7.937E+00 7.841E+00 7.787E+00 7.778E+00 
IY= 13 9.513E+00 9.524E+00 9.161E+00 8.915E+00 8.788E+00 
IY= 12 1.207E+01 1.208E+01 1.090E+01 1.012E+01 9.709E+00 
Iz@= 11 1.256E+01 1.254E+01 1.113E+01 1.022E+01 9.738E+00 
TY= 10 I 3256701 1.327E+01 1.147E+01 1.034E+01 9.750E+00 
l1Yy= 9 1.464E+01 1.452E+01 Pakese.e) 1.033E+01 9.563E+00 
Iy= 8 1.447E+01 1.435E+01 1.179E+01 1.022E+01 9.402E+00 
ie 7 1.439E+01 1.426E+01 1.174E+01 1.010E+01 9.233E+00 
IY= 6 1.298E+01 ez 76Et+O! 9.974E+00 8.171E+00 7.254E+00 
Iy= 5 1.232E+01 1.209E+01 9.194E+00 7.466E+00 6.610E+00 
Iy= 4 7.458E+00 7.186E+00 4.339E+00 3.453E+00 3.171E+00 
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IY= 3 2.818E-05 
IY= 2 -9.134E-01 
IY= 1 -*3.167E-03 
IX= 1 
Iv= 15 5 .653E+00 
IY= 14 7.782E+00 
IY= 13 8.656E+00 
IY= 12 9.296E+00 
IY= 11 9,231E+00 
IY= 10 9.113E+00 
Iv= 9 8.751E+00 
IY= 8 8.533E+00 
i= = 8.280E+00 
IY= 6 6.266E+00 
IvyY= 5 5 .671E+00 
IYy= 4 2.876E+00 
IY= 3 1.243E+00 
IY= 2 ones 
IYy= 1 .235E-01 
IX= oe 
FIELD VALUES OF W1 
ITY= 16 -2.942E+00 
IY= te -6.774E-Ol 
IY= 14 =. SLOE+00 
IY= 13 4.553E+00 
IY= 12 6.727E+00 
IY= 11 6.955E+00 
IY= 10 7.133E+00 
IY= 9 7 .025E+00 
IY= 8 §.555E+00 
hy= 7 4 .606E+00 
IY= 6 6.964E-05 
IY= 5 2.638E-04 
Ty= 4 2.840E-04 
Iv= 3 8.589E+01 
IY= 2 8.535E+01 
Ly= 1 8.678E+01 
IX= a 
IY= 16 -2.318E+00 
IY= 15 -6.671E-01 
IY= 14 Zee 0 
IY= 13 6.008E+00 
IY= 12 1.051E+01 
IY= 11 1.136E+01 
IY= 10 1.302E+01 
IY= 9 i; SO3E+O1 
IvyY= 8 1.614E+01 
IvY= 7 Pees 7Et+Ol 
LyY= %% 2.451E+01 
IY= 5 2.843E+01 
IyY= 4 3.479E+01 
Ey= 3 4.035E+01 
tY=p) 2 4.356E+01 
Iy=) 1 4.379E+01 
IXx= 6 
FIELD VALUES OF KE 
IY= 16 1.973E+02 
IY= 15 2.567E+02 
Ivy= 14 5.500E+02 
TY= 13 1.046E+03 
IY= 12 1.767E+03 
Iv= 11 1.910E+03 
IY= 10 2.209E+03 
IY= 9 2.615E+03 
IYy= 8 2.832E+03 
IY= 7 3.444E+03 
IY= 6 4.522E+03 
TEz= 5 5.295E+03 
IvyY= 4 6.637E+03 
IY= 3 7.055E+03 
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-841E+01 
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2-419E-08 
9.804E-01 
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“2, SOU Ea oO 
-6. Ooemaen 


2.22Z2E+00 
5 .500E+00 
9.514E+00 
1.024E+01 
1.166E+0O1 
1.344E+01 
1.442E+01 
1. (3ZE or 
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2 Gooey 


4. -968E+01 
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.851E+02 
. 5aZE+03 
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. Sazerus 
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~995E+02 
~547E+02 
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.874E+03 
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257 SE+03 
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WWWWNNNRRRrWOON-: 
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~152E+03 
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-870E+03 
~-128E+03 
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Sens 


654E+03 


IX= 6 
BEEED VALUES OF EP 


Iv tt 
LX= 
IY= 16 
Iy= 15 
IY= 14 
me 13 
IY= 12 
ty= 11 
TY= 10 
i= 9 
a 6 
= 7 
IY= 6 
lee 5 
IvY= 4 
IY= 3 
iz 2 
ma 1 
IY= 16 
IY= 15 
TY= 14 
m= 13 
He 12 
y= 11 
IY= 10 
i= 9 
IYy= 8 
IY= 7 
IY= 6 
Iy= 5 
IY= 4 
—- 3 
Ba- 2 
re | 
TX= 
IY= 16 
me 15 
IY= 14 
m= 13 
Be 12 
I 11 
TY= 10 


_ 
KG 
i 

HN WB O10) ~1 00 0 


MDM MINRPRRRY O~AIDWeE BW OMFRDWN BNR RPRrRrOWFE Pl 
- 


~-192E+02 
.939E+02 
~474E+03 
~910E+03 
-O08E+03 
-0O06E+04 
»2597E+04 
.654E+04 
.876E+04 
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ISWP EP Hl 
S.100E+02 3.604E+14 1.418E+13 
5.250EtO02 3.470Et14 2.645E+13 
5.400Et02 3.406E+14 2.557E+13 
VARIABLE Pl Ul Vi1 Wl KE 

MINVAL= 1.831E+01 2.193E+01 2.207E+01 2.318E+01 2.832E+01 
MAXVAL= 1.853E+O1l 2.214E+01 2.227E+01 2.359E+01 2.865E+01 
VARIABLE EP Hl 
MINVAL= 3.346E+O1l 3.028E+01 
MAAZVAL= 3.352E+01 3.091E+01 
1.0088. ...t+ec cet. s eee ot, ot, ot. eolenen esa V 
0.90 + H 
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ME 9-lOE+O2 MAX= 5.40E+02 
KAKKAAKAKKKKKKKAKAKKKAKAKKARKKARKKAK 
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eS 
OmM++++++4+'0 
OM+++++A++ 


Ts ABSCTSSA IS ope 
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A-t OV + 
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DATA FOR RE-STARTS AND PLOTTING SAVED ON DFO9 


RUN NO. 1 ENDED AT ISWEEP= 990 AND ISTEP= i 
RRAKARKKAAKKAKKAAKKAKKAKKKKKKAKKAKAKKAKAKKKAKKAKKAKAKKKKAKAKKKKK AR 
SATLIT RUN NUMBER = 1 

MACHINE-CLOCK TIME OF RUN = ees eCOND a. 


TIME/ (VARIABLES*CELLS*TSTEPS*SWEEPS*ITS) = 0.000E+00 
KRARRAARKAKAAKRKKAAKAKKRARAAARAKARKARKKAKRAKRKKRRKRKKKRARRERKARR 


NORMAL STOP REACHED IN PROGRAM 


ie 
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